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FOREWORD 


This  report  was  prepared  in  the  Applied  Mathematics  Research 
Branch  of  the  Aeronautical  Research  Laboratory,  Air  Force  Research 
Division,  Air  Research  and  Development  Command,  Wright- Patter  son 
Air  Force  Base,  Ohio.  The  problem  stems  from  a  proposed  electronic 
method  for  extracting  information  from  very  low-contrast  photographs, 
through  grain  counts  of  the  enlargements,  and  originated  in  the  Solid 
State  Physics  Research  Branch  of  the  Aeronautical  Research  Labora¬ 
tory.  The  authors  gratefully  acknowledge  numerous  discussions  with 
Dr.  Lee  Devol  and  Mr.  Radames  Gebel  on  the  basic  physical  question 
which  has  been  presented  in  WADC  TN  58-110,  "Electronic  Contrast 
Selector  and  Grain  Spacing  to  Light  Intensity  Translator  for  Photo¬ 
graphic  Enlargements",  by  R.  Gebel.  The  work  presented  herein  is  an 
attempt  at  a  suitable  mathematical  formulation  as  a  basis  for  proper 
analysis  of  such  grain  counts.  The  authors  also  acknowledge  the  assist¬ 
ance  of  Mrs.  Martha  Elmore  and  Mr.  James  Caslin  in  the  preparation 
of  some  of  the  tables  and  graphs.  Appendix  III  on  computational  procedures 
developed  for  formulas  given  in  Section  7  (to  be  used  with  the  Burroughs 
E-101  Electrodata  Electronic  Computer)  was  written  by  Mr.  Caslin.  The 
work  was  carried  out  under  Project  7071,  Task  70437,  "Methods  of 
Mathematical  Physics". 
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ABSTRACT 


The  accuracy  of  the  evaluation  of  a  photographic  plate  is  limited 
by  its  grain  structure.  One  approximates  the  value  for  the  light  density 
at  a  given  point  by  the  average  light  density  in  a  small  area  (the  "test 
are")  surrounding  the  point.  This  paper  establishes  confidence  limits 
for  evaluation  procedures  of  this  kind.  It  is  assumed  that  the  grains  on 
the  photographic  plate  arise  in  independent  random  processes  controlled 
by  the  local  density  of  the  light  flux.  In  the  evaluation  procedure  one 
counts  the  number  of  grains  in  the  test  area.  Generalizing  the  method 
one  attaches  a  weight  factor  to  each  grain  depending  upon  the  grain 
position  within  the  test  area  and  then  determines  the  sum  of  the  weight 
factors  for  the  grains  found  in  the  test  area.  By  such  a  procedure  one 
can  determine  quantities  related  to  the  light  density,  e.  g.  the  density 
gradient;  one  can  also  scan  for  patterns  of  a  special  kind,  e.  g.  a  sudden 
jump  of  the  light  density.  For  measurements  of  this  kind  probability 
theory  predicts  the  expected  value  and  the  variance  in  terms  of  the  light 
density  and  the  chosen  weight  function.  There  are  two  kinds  of  errors 
in  the  measurement  process:  errors  due  to  the  non-vanishing  size  of 
the  test  area,  and  errors  due  to  the  randomness  inherent  in  the  process 
of  grain  generation.  The  variance  due  to  errors  of  both  kinds  must  be 
minimized.  The  treatment  of  these  questions  is  shown  in  a  number  of 
examples  of  increasing  complexity.  Moreover,  this  report  investigates 
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the  probability  of  making  an  error  if  one  tries  to  discriminate  between 


two  known  light  densities  on  the  basis  of  grain  counts,  and  it  also  examines 
a  step-wise  method  for  carrying  out  such  counts. 
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INTRODUCTION 


In  order  to  gain  maximum  information  from  a  photograph  one  may 
enlarge  it  to  the  extent  that  the  individual  grains  show  up.  The  density 
of  the  grains  is  proportional  to  the  flux  of  light  hitting  a  particular  area. 
The  grains  of  the  photographic  plate  arise  by  a  random  process;  there¬ 
fore,  the  statistical  variations  in  the  density  of  them  set  a  limit  to  the 
evaluation.  The  interpretation  of  the  photograph  will  then  amount  to 
somehow  forming  average  densities;  also  one  may  try  to  recognize 
patterns  or  some  main  features  of  the  picture.  This  process  is  usally 
done  by  inspection,  perhaps  with  some  additional  photographic  techn:.ques, 
such  as  high  contrast  prints.  The  question  can  be  asked  whether  this 
process  can  be  done  by  a  photoelectric  scanning  process.  This  may  be 
desirable  for  several  reasons.  The  photoelectric  processes  might  reveal 
information  which  is  hard  to  recognize  with  the  eye.  Whether  this  is 
true  dep  nds  upon  the  performance  of  the  eye,  and  this  question  does  not 
concern  us  in  this  report.  The  photoelectric  evaluation  would  not  require 
human  judgment,  and  thus  it  is  not  subject  to  fatigue  and  also  it  can  be 
done  more  quickly.  It  is  felt  that  the  second  point  is  sufficiently  important 
to  justify  a  closer  study  of  this  possibility. 

We  ask  in  this  report  which  conclusion  about  the  energy  density  of  the 
electromagnetic  waves  can  be  drawn  from  the  distribution  of  grains  on 
given  photographic  plates.  Local  measurements  are  impossible  because 
of  the  grain  structure  of  the  plate,  therefore  one  will  try  to  characterize 
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the  desired  light  density  at  a  given  point  by  an  average  -  possibly  a 
weighted  average  -  of  the  grain  density  over  an  area  surrounding  this 
point.  This  paper  treats  this  process  from  a  probabilistic  point  of 
view.  Confidence  limits  are  established  for  different  measuring  pro¬ 
cedures  and  it  is  shown  how  to  choose  a  weight  factor  characterizing 
the  measurements  in  such  a  manner  that  the  error  limits  are  minimized. 

By  a  suitable  choice  of  the  weight  factor  it  is  possible  to  determine, 
beside  the  light  density, other  quantities  related  to  it,  e.  g.  its  gradient. 
Furthermore, it  is  possible  to  search  for  special  patterns,  e.  g.  for  a 
jump  of  the  light  density,  or  for  narrow  lines  of  a  higher  density. 

From  a  technical  point  of  view  the  question  of  discriminating  between 
two  fixed  light  densities  may  be  of  interest.  Here  one  will  define  a  certain 
cut-off  count  for  the  number  of  grains,  and  ascribe  to  all  areas  with  a 
count  below  the  cut-off  the  lower  light  density,  to  all  others  the  higher 
one.  For  this  procedure,  formulae  for  the  probability  of  error  are  giver. 

From  a  practical  point  of  view,  one  might  be  inclined  to  determine 
first  for  subareas  of  the  test  area,  whether  they  belong  to  the  higher  or 
to  the  lower  light  density.  The  majority  of  assignment  of  one  kind  or  the 
other  in  the  test  area  will  then  determine,  whether  the  test  area  belongs 
to  the  higher  or  to  the  lower  count.  The  probability  of  error  for  this 
procedure  and  for  a  direct  count  for  the  entire  test  area  is  also  investi¬ 
gated. 
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1.  Model  Connecting  the  Grains  of  the  Photographic  Plate  With  the 


Light  Flux 

For  our  theoretical  analysis  we  need  a  model  for  the  interaction 
of  light  with  the  photographic  plate.  The  light  manifests  itself  in 
photons,  i.  e.  at  discrete  points  of  the  photographic  plate.  Not  every 
photon  gives  a  chemical  reaction  on  the  plate;  a  major  portion  of  the 
light  passes  through  without  indication.  But  some  photons  modify 
certain  molecules.  The  grains  which  we  see  are  clusters  of  changed 
molecules  which  grow  because  of  the  development  process  around  the 
molecules  of  the  photographic  emulsion  that  have  been  changed  by 
the  photons.  The  subsequent  computations  are  based  on  the  follow¬ 
ing  model  for  the  creation  of  grains:  For  a  first  example  assume 
the  energy  density  of  the  light  waves  to  be  constant.  Let  us  assume 
that  we  can  stretch  the  time  scale  in  such  a  manner  that  one  grain 
arises  after  the  other.  Then  the  grain  arising  first  will  be  found  at 
some  spot  of  the  plate  and  all  locations  are  equally  probable.  The 
probability  of  finding  a  grain  in  a  given  area  is  equal  to  this  area 
divided  by  the  area  of  the  whole  plate.  The  second  grain  arises  in 
exactly  the  same  manner  and  its  location  is  not  influenced  by  the 
location  of  the  first  grain.  The  sa'ne  holds  for  all  grains  that  arise 
subsequently.  One  says  that  the  arising  of  grains  represent  inde¬ 
pendent  events. 

If  the  light  density  varies  across  the  plate,  then  the  probability 
of  finding  the  first  grain  in  a  given  area  is  equal  to  the  integral  of 


3 


the  light  flux  for  the  area  considered,  divided  by  the  integral  for 
the  light  flux  over  the  whole  plate.  Again  the  production  of  any 
one  grain  is  assumed  to  be  independent  of  the  production  of  all 
other  grains.  Furthermore  it  is  assumed  that  all  grains  have  the 
same  darkness  and  can  be  considered  as  points,  i.  e.  if  one  counts 
grains  in  a  given  area,  a  grain  is  either  inside  or  outside  of  that 
area. 

2.  Some  Results  of  Probability  Theory 

For  our  future  discussion  some  results  of  probability  theory 
will  be  needed.  They  are  derived  in  Appendix  I.  Here  they  will 
be  quoted  and  described  in  a  self-contained  manner. 

The  number  of  photons  that  act  on  the  photographic  plate  and  the 
number  of  grains  produced  are  considered  as  very  large.  For  the 
measurement,  a  small  area  called  "test  area"  of  the  photographic 
plate  is  considered  and  from  the  number  and  the  arrangement  of 
the  grains  certain  conclusions  on  the  density  of  the  incoming  light 
are  drawn. 

Let  us  introduce  on  the  phographic  plate  a  Cartesian  system  of 
coordinates.  The  probability  that  a  grain  is  produced  in  a  certain 
area  is  considered  to  be  proportional  to  the  light  flux  across  the 
area  multiplied  by  the  exposure  time.  Let  the  density  of  the  light 
flux  (the  energy  of  the  electromagnetic  waves  per  unit  area)  be 
given  by  p(x,y).  If  one  makes  a  great  number  of  photographs  of 
the  same  object  and  counts  in  all  photographs  the  number  of  grains 


in  the  same  test  area  A,  then  the  counts  will  in  general  not  be  the 
same.  Let  N  be  the  conceptual  average  formed  from  an  infinitely  great 
number  of  such  counts.  (The  number  N  is  also  denoted  as  "  expected 
value" )•  A  factor  of  proportionality  which  contains  the  exposure  time 
and  the  sensitivity  of  the  plate  may  be  denoted  by  k  such  that 

N  =  k  //  p  (x,  y)  dxdy  (1) 

A 

is  the  average  number  of  grains  to  be  expected  in  A. 

Let  0(x,  y)  be  a  weighting  function.  The  choice  of  this  weight¬ 
ing  function  determines  the  measurement  that  is  carried  out.  The 
individual  grains  found  in  the  test  area  in  a  specific  experiment  may 
be  numbered  as  1, 2,  ...  i  . . .  n  (n  will  be  close  to  N  but  need  not  be 
equal  to  N);  the  coordinates  of  the  individual  grains  are  (x^.y^)  .  . . 
^Xi’yi^  ”  *  ^he  *unct*on  0  (x.  y)  then  yields  the  weight  0(x^,y.)  for 
grain  i  in  the  position  (x.,y.) .  A  measurement  then  determines  the 
expression 


Y=Z0(xiiy.)  (2) 

where  the  summation  is  to  be  extended  over  all  points  that  lie  in  the 
test  area.  For  a  fixed  area  A  and  a  given  photographic  plate,  Y  will 
have  a  fixed  value.  If  one  evaluates  a  great  number  of  plates  which 
have  arisen  from  the  same  experiment,  i.  e.  from  the  same  function 
p(x,y)  with  the  same  constant  k,  one  will  obtain  different  values  of 
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Y.  Let  p(u)  be  the  probability  density  of  Y  taking  on  values  u.  Fig.  1 
shows  such  a  probability  density  curve.  The  area  under  the  curve  is 
the  probability  of  finding  a  test  result  between  u  =  •  «•  and  u  =  + «°.  This 
probability  is  1.  The  average  value  of  Y  is  defined  by 


+«° 

my  =  /  u  P<u)  du 


(3) 


In  our  measurements  the  quantity  that  is  to  be  determined  will  be 
approximately  proportional  to  the  average  value;  but  naturally,  indi¬ 
vidual  measurements  will  give  results  that  are  scattered  around  the 
average  value.  To  characterize  the  width  of  the  scatter  one  introduces 
the  variance  which  is  defined  by 


+«0 


Vy  =  /  (u-My  )  p(u) du 


(4) 


To  show  how  the  variance  is  connected  with  the  width  of  the  curve, 
let  us  first  assume  that  the  probability  is  constant  over  values  Y 
which  deviate  from  the  average  by  at  most  b.  (In  other  words  b  is  the 
half  width  of  a  rectangular  distribution).  Then  the  probability  density 
is  l/2b  and  one  obtains  for  the  variance 


+b  2 

V  =  —  /  (u  -  M  r  d(u  -  M  )=-^— 
Y  2b  -b  Y'  'u  "V  3 
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i.  e.  in  this  case 


b  =  r  ^  (5) 

where  r  =  */3~ 

Frequently,  e.  g.  if  the  number  of  points  in  the  test  area  is  large, 
the  distribution  will  be  close  to  normal,  as  can  be  shown  by  means  of 
the  Central  Limit  Theorem  (see  [1]  page  180).  Then,  since  the  distri¬ 
bution  extends  to  infinity,  b  can  no  longer  be  defined  as  the  half  width 
of  the  area  which  contains  all  the  points.  It  is  then  considered  as  the 
half  width  of  an  area  that  contains  a  fraction  (1  -  a)  of  all  points.  We 
can  use  (5)  again;  r  is  given  as  a  function  of  a  in  the  following  table 


a 

1  -  a 

T 

.317 

.683 

1 

.045 

.955 

2 

.003 

.997 

3 

.  100 

.900 

1.64 

.050 

.950 

1.96 

.010 

.990 

2.58 

The  value  of  t  which  we 

obtained  previously  for  a  rectangular 

distribution  then  corresponds  to  an  a 

of  about  0.  1,  i.  e.  90%  of  the 

measurements  will  fall  within  the  region 

1.6  n/Vy 
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It  is  not  impossible  to  determine  the  probability  density  curves 
for  the  experiments  to  be  described  in  this  report.  They  may  be 
obtained  from  the  characteristic  functions  given  in  Appendix  I.  Here 
we  mention  only  two  results.  The  average  value  of  infinitely  many 
measurements  (i.  e.  for  an  identical  test  area  in  infinitely  many  photo¬ 
graphic  plates)  is  given  by 

My  =  k //  p(x,y)  0  (x,y)dxdy  (6) 

A 

The  variance  is  given  by 

V  =  k  //  p  (x,  y)  02  (x,  y)  dxdy  (?) 

A 

In  making  measurements  one  will  define  6  in  such  a  manner, 
that  the  average  value  My  is  directly  related  to  the  quantity  which 
one  wants  to  find.  Then  one  can  compute  from  (7)  the  variance 
incurred  in  such  a  measurement  and  determine  from  it  probable 
bounds  for  the  errors  which  may  arise,  if  one  uses  Y  as  an  approxi¬ 
mation  for  My.  Examples  will  be  given  in  the  next  section. 

Later  another  result  of  probability  theory  will  be  needed.  The 

quantity  as  defined  by  the  right  hand  side  of  (6)  is  in  general  only  an 

approximation  to  the  quantity  which  is  to  be  determined.  One  reason 

is,  that  the  form  of  (6)  is  an  average  over  an  area,  while  we  are 

frequently  interested  in  local  values.  To  fix  the  ideas  take  the  local 
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value  of  p  (x,  y).  This  local  value  of  p  must  be  approximated  by  an 
integral  My  over  an  area.  In  general  this  average  My  will  not  be 
equal  to  the  value  of  p  at  the  desired  point.  The  difference  depends 
upon  the  character  of  the  functions  p  and  0,  and  it  will  vary  with 
x  and  y;  basically  it  depends  upon  the  properties  of  the  photograph. 
Thus  the  quantity  My  -  even  if  suitably  defined  -  is  only  an  approxi¬ 
mation  to  p  and  the  difference  My  -  p  possesses  a  variance,  which 
we  shall  denote  by  .  In  estimating  My  by  counting  the  number  of 
grains  within  the  area  and  calculating  Y  another  error  is  incurred 
which  has  a  variance  V  .  Now  probability  theory  gives  the  following 
result:  If  the  two  probabilistic  processes  are  independent  then  the 
result  will  have  a  variance  given  by  the  sum  of  the  two  variances. 

3.  Examples 

Some  examples  may  show  how  these  results  are  applied  to 
specific  situations,  (a)  Let  us  assume  that  p  is  constant  (but  unknown) 
and  that  we  want  to  determine  its  value.  Taking  a  rectangular  test  area 
with  sides  2a  and  2b  one  finds 

+a  +b 

my  =  kP  /  /  6  <x*  y)  dxdy 

-a-b 


a  b 

k  f  f  6  (x,  y)  dxdy 
-a  -b 
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The  last  equation  shows  the  relation  between  the  desired  value  p 
and  the  average  value  p  is  proportional  to  and  the  factor  of 

proportionality  is  known.  The  variance  is  given  by 

+a  +b  2 

Vv  =  kp  /  /  0  (x,  y)  dxdy 

Y  -a  -b 

In  these  formulae  6  is  still  arbitrary.  Let 

0  =  1 

then  My  =  kp4ab 

My  is  then  the  average  number  of  points  N  in  the  area.  The  variance 
is 

Vy  =  kp4ab  =  My  =  N 

Using  Eq.  (5)  for  the  width  of  the  probability  curve  one  finds  as 
measure  for  the  relative  error 

e  =  JL  =  -x_  (8) 

>JfT 

Here  t  must  be  taken  from  the  table  given  in  Section  2.  If  a  = .  003 
and  N  =  10,  000  one  would  obtain  a  relative  error  in  the  order  of  3%; 
for  N  =  100  one  would  obtain  30%  relative  error.  For  a  =  .  10  the 
relative  error  with  1000  points  is  5.2%.  This  may  be  considered  as 
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tolerable.  About  1000  points  in  the  test  area  appears  to  be  the  minimum 
that  should  be  used. 

One  might  try  to  reduce  the  variance  by  choosing  another  function  6. 
Unfortunately  6  =  constant  is  the  best  choice.  Let  us  pose  the  problem 
of  minimizing  the  variance  while  the  average  value  is  kept  constant: 

//  0(x,y)  dxdy  =  constant 
A 


//  0^  (x,  y)  dxdy  =  minimum 
A 

This  is  a  simple  problem  of  the  calculus  of  variations.  Introducing 
a  Lagrangian  multiplier  X  cne  obtains  an  equivalent  problem: 

2 

//  {0  (x,y)+ X0(x,y)}  dxdy  =  minimum 
A 

Hence 


//  6  0  {  2  0(x,  y)+  X  }  dxdy  =  0 
A 

or 


0  =  - 


X 

2 


constant 


In  computing  the  relative  error  the  choice  of  the  constant  in  the 
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last  equation  is  unessential.  We  set  0  =  1,  and  thus  find  that  the  above 
choice  of  0  was  optimal. 

(b)  Let 

p  =  p0[1  +  al  I  +bl  b  1  1  (9) 

we  try  to  determine  p  at  x  =  0,  y  =  0.  Choosing  0  =  1  we  obtain 

MY  =  kpo  4ab 

Thus  is  proportional  to  at  the  origin,  which  shows  that 
this  choice  of  0  is  suitable  for  the  desired  measurement.  The 
variance  is  found  to  be 

Vkpoi.{<1  +  al  I  +b1i)dxdy 
=  4*P0ab 

Again  these  quantities  can  be  expressed  by  means  of  the  average 
number  of  points  N  in  the  test  area,  and  one  obtains  the  same  result 
for  the  relative  error  as  in  (8. ) 

Let  us  again  consider  the  minimum  problem  for  9 .  The  function 
0  (x,  y)  must  be  chosen  such  that  only  the  constant  term  in  the  above 
expression  for  p  will  give  a  contribution  to  •  Let 
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Pj  (x,  y)  =  1 


P2(x,y)  =  x/a 

P3(x.y)  =  y/b 


Then  we  must  determine  0  (x,  y)  in  such  a  manner  that 


+a  +b 

/  /  9(x,  y)  p.  (x,  y)  dxdy  =  1  (10a) 

-a  -b 


V)  0(x,y)p  (x,y)dxdy  =  0  (10b) 

-a  -b  L 

+a  +b 

/  ^  0  (x,  y)  p3  (x,  y)  dxdy  =  0  (10c) 

+a  +b  2 

-J  pi +  ai  p2 +  bi  P3^9  (x*  y)dxdy=  minimum  (lOd) 


By  variational  calculus  one  obtains 


0(x.y)  = 


X1  P1  *  X2  p2  *  Sp3 

pl+al  p2+ V3 


where  and  \  must  be  determined  from  conditions  (10a, b,c). 

By  chance  this  problem  possesses  a  simple  solution.  If  we  set 


13 


one  finds  0  =  1.  The  conditions  (10b)  and  (10c)  are  automatically 
fulfilled.  The  condition  (10a)  is  unessential,  it  could  be  fulfilled  by 
choosing  a  suitable  constant  for  6,  instead  of  6  =  1.  But  such  a 
constant  would  drop  out  in  the  determination  of  the  relative  accuracy  e  . 

Thus  we  find,  that  even  for  a  p  surface  that  possesses  a  gradient, 

9  =  1  is  the  best  choice  for  the  weight  function. 

(c)  Going  one  step  further  one  might  consider  a  function  p  which 
is  represented  in  a  certain  area  by 


P=  P0  (  1+ajX  +  a2y  +  a3x 


a4xy 


+  a, 


V2) 


To  obtain  sumple  results  we  assume  that  p  is  nearly  constant.  Then 
by  calculus  of  variations  one  finds 


0 


~  X,,  +  x  +  \,  y  + 
12  3  4 


2 


x 


+ 


xy  +X6 


2 

y 
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with  the  additional  conditions 


//0x2dxdy  =  O  (lid) 

A 

/  /  0xy  dxdy  =  0  (lie) 

/  /  0  y2  dxdy  =  0  (Ilf) 

For  reasons  of  symmetry  one  finds  from  (lib)  =  0,  from 
(11c)  X.3=0,  from  (lie)  X.  0.  If  we  consider  a  region  A  which 
does  not  change  if  the  x  and  y  axis  are  interchanged,  one  further¬ 
more  finds  X.  .  =  \6  .  Let  X  =  1.  We  now  have 

0  =  l  +  X6(x2  +  y2) 
with  the  additional  conditions 

/  /  0  x2  dxdy  =  0 
A 

/  /  0  y2  dxdy  =  0 
A 

or  /  /  0  (  x2  +  y2  )  dxdy  =  0 

A 


/  /  0  dxdy  =  1  (11a) 

A 


/  /  0xdxdy  =  0  ( lib) 

A 


/  /  0  y  dxdy  =  0  (11c) 

A 
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If  we  consider  a  circular  region  of  radius  1,*  then  this  con¬ 


dition  will  give 


1  2  3 

/  (  1  +  X ,  r  )  r  dr  =  0 

0  6 

Xa=-3/2 


Thus 


Q  ,  3  2 

0=1-Tr 


With  this  value  one  obtains  for  the  average  measurement 


My  =2irkP()  /  (  1  -  |  r2  )rdr  =  |  TrkP() 


The  variance  is 


1  3  2  2  1 
Vy  =  2irkp0  /  (  1  -  j  r  )  rdr  =  -  irkp( 


The  number  of  points  in  the  region  is  N  =  irkpg 
Thus  the  average  of  the  measured  quantity  is  N/4 
and  the  variance  is  N/4 
The  relative  error  is 


€  - 


2t 

n/N 


♦Since  the  scale  of  the  x  and  y  coordinates  should  be  unessential,  it 

can  always  be  chosen  in  such  a  manner,  that  at  the  boundaries  x=  1,  or 
y  =  1,  or  r  =  1  (whichever  is  convenient). 
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If  one  would  have  omitted  the  conditions  that  the  result  has  to 

2  2 

be  sensitive  against  contributions  of  x  and  y  in  p  then  the 
relative  error  would  have  been  1/2  of  this  amount(  viz 


€  = 


For  comparison  the  same  computation  may  be  carried  out  for 
a  square  region. 

.  ,  .2  2  . 

0  =  l+\6(x  +  y  ) 

X. .  is  determined  from  the  condition 


Thus 


Thus 


/  /  0  x  dxdy  =  0 
A 


/  /  [  x2  +\,  (  +  x2  y2  )  ]  dxdy  =  0 

A  6 

6  i4 


.  ,  15  ,  2  ,2  . 

0=1-  ~  (x  +y  ) 

11  11 
M  =4kp  //  0dxdy  =  4kp  //  [  1  - 
0  0  0  0 


15  .  2  2..  8  , 

—  (x  +y  )]dxdy  =  j  kpQ 


The  number  of  points  in  this  region  is  4kp^  =N 
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Thus 


The  variance  is  found  to  be 

1  1  15  2  2  2 

VY  =  4kp0  q  q  ^  1  ‘  14  (  X  +  Y  ^  dxdy  =  7  kp0  =  MY 
Here  the  relative  error  is  given  by 

VT5t 

i.  e.  for  the  square  the  relative  error  is  somewhat  smaller  than  for 

a  circle  if  the  same  number  of  points  is  used. 

The  limitations  of  this  procedure  are  quite  obvious.  While  in 

principle  it  is  possible  to  exclude  by  constraints  the  influence  of 

additional  terms  in  a  development  for  p  each  constraint  may  make 

the  function  0  more  eneven  and  thus  increase  the  variance.  We 

* 

shall  take  up  this  question  later. 

To  show  another  variation  of  this  technique  let  us  try  to  determine 
a  procedure  for  finding  the  first  derivative  of  a  given  function  in  the  x 
direction. 

Given 


P  =  PQ(  1+aix  +  a2y) 


(12) 
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Assume  again  that  a^  and  a^  are  small,  so  that  p~  p^.  Then 
an  optimal  expression  6  has  the  form 

0=  Xj  +  X^x  +  X^y 
subject  to  the  conditions 

/  /  0  x  dxdy  =  1 

A 

/  /  0  dxdy  =  0 

A 

/  /  0  y  dxdy  =  0 

A 

The  second  condition  gives  immediately  X^=  0,  the  third  one 
X3=  0.  Let  the  test  area  be  a  square  of  side  2a.  The  first  con¬ 
dition  yields 

0=  ^-2-  x  (13) 

4a4 


One  then  finds 


+a  +a 

M  =  /  /  0(x)  p  (x)  dxdy  =  a^pg 
-a  -a 

Let  N'  =  kpg  be  the  expected  number  of  points  per  unit  of  area 
(cf.  (1)).  Then 

My  =  ajN* 
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Hence 


a 


1 


(14) 


The  variance  is 


+*+?  2  3kp0  3N' 

VY  =  /  /  0  (  x,  y)  p  (  x,  y)  dxdy  =  7  —  =  4  “4 

Y  -a -a  a  a 


(15) 


According  to  the  remarks  in  Section  2  the  half  width  of  a  region 
which  contains  a  fraction  1  -  a  of  all  measurements  is  tn/V^  • 
The  uncertainty  in  the  measured  values  of  Y  can  be  expressed 


as  an  uncertainty  in  the  values  of  a^  .  If  the  measurement  Y 
deviates  from  by  WVY  then  one  obtains  from  (14)  as  devi¬ 


ation  in  a 


1 


Aal  = 


T'iyJ 

N' 


W3_ 

2aZ'^?^, 


For  a  =  0.  1  (cf.  Section  2)  one  finds  t  =  1.  6  and 


a  2.76 

Aai  =  7T- 


2a*N/NT 


(16) 


For  a=  1  the  average  number  of  points  in  the  test  area  is  4N1. 
Assume  4ir  =  900.  Then 


Aaj  =  0.  092 
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The  significance  of  can  be  seen  from  Eq  (12).  The  uncertainty 
in  a^  just  found  is  rather  large.  Eq  (16)  shows  how  the  size  of  the 
test  area  will  influence  the  accuracy. 

4.  More  general  distribution  of  p 

In  previous  considerations  we  assumed  a  certain  form  of  the 
functions  and  then  determined  the  accuracy  with  which  the  value  at 
a  given  point  could  be  obtained.  The  result  was  invariably  that  the 
region  whould  be  taken  as  large  as  possible.  A  limit  arises  naturally 
by  the  fact  that  the  analytic  expressions  for  p  used  here  will  repre¬ 
sent  the  actual  p  (x,y)  in  a  limited  region  only.  For  functions  p  as 
they  occur  in  practice  the  expressions 

MY  =  /  /  6(x,y)p(x.  y)dxdy 
A 

are  only  approximations  of  the  desired  quantities  p  or  grad  p. 

To  extend  the  techniques  used  so  far  to  more  and  more  complex 
analytical  expressions  for  p  would  be  rather  useless,  for  with  this 
process  the  variance  will  increase.  Therefore  we  must  accept  the 
fact  that  the  average  values  which  we  determine  are  only 
approximations  of  the  desired  quantities. 

In  the  following  discussions  we  shall  use  for  the  determination 
of  p  at  a  given  point  an  approximation  which  is  based  on  the  assump¬ 
tion  that  in  the  area  considered  p  is  a  linear  function  of  x  and  y. 
Correspondingly  0=1.  The  value  p  in  the  center  of  the  area  is  then 
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approximately  proportional  to 


+a  +a 

MY  =  /  /  p(x,y)dxdy 
-a  -a 

where  the  origin  of  the  xy  -  system  lies  at  the  center  of  the  region 
considered. 

On  a  given  photograph  there  exist  lines  of  constant  p  .  Let  us 

consider  along  one  of  these  lines,  say  p  =  p^,  a  great  number  of 

equidistant  points.  These  points  may  be  numbered  as  1.2...  j... 

their  coordinates  will  be  (x^,  y^ )  ...  (x.,  y. )  . . .  .  We  then  obtain 

J  1 

for  the  value  for  point  j 

y.+a  y.+a 

j  j 

M  =  /  /  p(x,y)dxdy 

’  ^  y.  -a  y .  -a 
J  J 

Since  the  value  of  will  vary  from  point  to  point  (although  all 
points  j  lie  on  a  line  p  =  constant),  one  may  regard  itself  as  a 
random  function.  For  this  random  function  there  exists  a  probability 
distribution,  which  gives  the  probability  density  for  to  take  on 
given  values.  For  our  measurements  these  values  would  be 
interpreted  as  proportional  to  Pq.  Obviously  this  probability  distri¬ 
bution  and  its  variance  is  determined  by  the  character  of  the  object. 
Intuitively  it  is  clear  that  for  an  object  for  which  p  varies  only  slowly, 
has  a  smaller  variance  than  for  an  object  with  rather  abrupt 
changes.  In  general  the  values  of  this  variance  must  be  estimated. 
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But  it  may  be  of  some  interest  to  see  how  this  variance  can  be 
determined  in  principle. 

Assume  for  this  purpose  that  for  the  vicinity  of  the  point 
(  x.,  y. )  described  above,  p  is  actually  of  the  form 


p  =  p0[l+a1  < 


Xj)  +  a2.j(y'yj)+a3,j( 


v2+ 


a4,j(y-V  J 


The  deviation  of  M  .  from  the  value  p.  is  easily  computed  as 

i ,  j  u 


MY,j  •  "(T  3  "0<a3  +  a4)a4 


The  quantity  a^  +  a^  can  be  considered  as  a  random  variable 

which  takes  on  the  values  a„  .  +  a .  .  for  the  points  considered  and 

3,  J  4,j 


possesses  a  variance  V.  We  assume  that  the  average  of 

(a3+a4) 


the 


random  variable  a„  +  a .  is 
3  4 


+  » 


M(a  +  a  )  =  !  (a3+a4)p(a3+a4)d(a3+  a4)  =  0 

J  *  "  ^ 


The  variance 


+  “ 


V(a,  taj"  !  (a3  +  a4>  p  (a3  +  a4>  d  (a3  +  a4» 
3  4  -  oo 


is  considered  as  known. 
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The  variance  of  about  the  correct  value  p^  is  then  given 
by 


V(M  ) 

y 


p(a3+a4)  d(a3+  a4) 


hence 


V(M  )  = 

y  9 


a 


(a3  +  a4) 


The  value  of  this  calculation  lies  in  the  fact  that  it  shows  how 
V(My)  can  be  computed  and  how  it  depends  upon  "a"  and  other 
known  quantities.  Naturally  for  other  assumed  representations  of 
p  ,  a  similar  computation  can  be  carried  out. 

V(My)  is  the  variance  of  the  average  measurement  from 
the  light  density  p^.  For  each  value  p^  of  the  light  density  occur- 
ing  in  a  given  picture  the  variance  V(M^)  can  be  determined. 

Since  4p^a^  is  the  number  of  points  N  in  the  region  considered, 
the  last  result  can  also  be  written  as 


N2a4 

V«My»=  — 


The  variance  which  occurs  in  the  measurement  of  p  has  been 
shown  in  Section  3,  Example  a,  to  be  N.  According  to  the  theorem 
quoted  at  the  end  of  Section  2  the  two  variances  must  be  added.  Thus 
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the  total  variance  of  the  measurement  of  Y  about  p.  is 


,t2  4 
N  a 

N  +  n  V.  .  , 

9  (a3  +  a4) 


The  measure  for  the  relative  error  is  then 


t|/n+  V.a  +a 

€  =  *  9  (a3+a^ 


Since  N  =  4p^  a  ,  one  has 


ll  2"  16  r~ 

rPoa,?P0  V<a,+ 


(a^+  a^)e 


4Pna 


To  obtain  best  conditions  this  quantity  must  be  minimized, 
i.  e.  one  must  determine  the  minimum  of 


Hence 


4p0  +  ii  V  a' 

- r  9  <va4> 


(a3  +  a4)p0 


This  equation  defines  the  best  size  of  the  region.  One  will 
notice  that  this  result  depends  upon  the  value  of  p^;  therefore, 
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•fc"  must  be  expected  to  vary  with  p^.  Among  the  values  of  "a" 
so  obtained  from  one  photograph  one  must  choose  a  suitable  value. 


Since  V.  .is  in  general  not  known,  this  discussion  is 

U3  V 

somewhat  academic,  but  it  shows  how  the  size  of  the  test  region 
is  connected  to  the  character  of  the  photograph. 

No  essential  changes  would  occur  if  instead  of  (12)  a  different 
form  of  p  would  be  assumed,  e.  g. 


p  =  p0(1  +  a1x  +  a2y  +  a3  |x|  +  a4  |y|  ) 


5.  Scanning  for  Special  Patterns 

The  examples  up  to  now  were  concerned  with  the  determination 
of  p  and  of  grad  p .  It  may  occur  that  certain  special  patterns  of 
the:  light  distribution  are  of  importance  and  one  might  try  to  detect 
them.  These  patterns  might  be  of  such  a  nature  that  they  would  be 
averaged  out  if  one  simply  determines  p  and  tries  to  recognize 
them  afterwards.  An  extreme  example  may  be  diffraction  rings  in 
a  telescope  due  to  the  finiteness  of  the  aperture.  They  may  be  rather 
completely  hidden  in  the  background  light.  More  practical  examples 
are  sharp  lines  like  streets  or  rivers  which  are  so  narrow  that  they 
become  rather  indistinct.  Also  the  question  may  arise  whether  along 
a  certain  line  the  density  jumps,  or  whether  there  is  a  gradient.  Such 
cases  are  covered  by  our  previous  investigations.  If  the  function 
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describes  the  pattern  which  is  to  be  recognized,  and  p.  and  p 
are  functions  which  should  be  suppressed,  and  if  p  is  nearly 
constant,  then  the  function  8  (x,  y)  has  the  form 

0  =  X  1  PX  +X2I>2  +X3P3 

and  the  X's  are  to  be  determined  from  the  conditions  that  the 
influence  of  p^  and  p^  must  be  suppressed. 

Assume  e.  g.  that  p^  is  the  distribution  due  to  the  Newton  ring, 
and  that  p^  =  1  is  the  background  light. 

Let 


P  =P0(1  +  aiPl) 


then  one  has  0  =  p  ^  +  X  p^ 


and  from  /  /  0  p  dxdy  =  0 
A 


one  obtains,  setting  p^  =  1 


x  =  ■  7  /  /  p  l  dxdy 

A  1 


i.  e.  0  is  equal  to  p^  minus  a  constant  which  makes  the  average  of 
0  over  A  zero. 

The  average  of  the  measurements  is  then 


MY=kpoai//(pi+  X  p2)  px  dxdy 
A 


(19) 
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The  variance  of  the  measurement  is 


VY  =  kp0  ^<pl+  X‘p2^2  (2°) 

A 

To  get  a  feel  for  the  orders  of  magnitude  let  us  assume  that 
the  area  considered  is  A,  and  that  the  function  p^  is  1  in  the 
area  Yq  A  and  is  zero  in  the  rest  of  the  region,  where  Yq  Is  a 
constant  between  0  and  1  (See  Fig  2.)  Further,  for  simplicity, 
let  a^  be  small.  Then 

//  p  dxdy  =  y0  a 
A 

and 

8  =  •»  1  -  Y0  p2 

The  average  measurement  will  give 

MY  =  kp0ai  //<pl’y0p2,pldxdy 
A 

=  kpQa1  Ay0(1-Yq)=  Naj  Yq(1"Yq) 

Let  N'  =  NYq  be  the  number  of  points  in  the  area  where  p^  ^  0. 
Then  My  =  a^l-  y0)  N' 
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The  variance  is 


VY  =kp0  Y0(1-V0)  =  NY0(1-Y0)  =  n'(1-Y0) 

We  proceed  by  analogy  to  the  treatment  of  the  density  gradient. 
The  half  width  of  a  region  in  which  a  fraction  (1  -  a  )  of  all  measure¬ 
ments  Y  are  contained  is  given  by  t  N/Vy  . 

One  may  ask  for  which  value  of  a^  the  half  width  equals  the 
expected  value  My  ;  one  finds 


(21) 


The  aim  should  be  to  make  this  value  small.  One  limited  possi¬ 
bility  is  the  choice  of  Yq  .  However,  if  we  keep  N'  fixed  and 
choose  Yq  small,  then  we  increase  the  total  test  area,  while  the 
area  belonging  to  the  pattern  is  the  same.  In  this  process  we 
would  encounter  additional  errors  due  to  the  nonuniformity  of  the 
background  field.  A  value  of  Yq  =  1/2  appears  reasonable.  Using 
T  =  1.  6  (corresponding  to  a  =  0.  1)  one  obtains 

1.  6 

ai  "  \fT7l  n/^T' 

For  a  =  0  1,  i.  e.  for  a  pattern  which  is  10%  darker  than  the 

2 

background,  the  pattern  must  contain  N'  =  =512  points  to 

w 
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meet  this  condition. 


This  approach  applies  to  patterns  of  any  kind.  One  might,  e.  g. 
try  to  recognize  a  straight  line  of  points  arranged  over  a  certain 
width.  In  this  case  one  would  consider  a  test  area  of  about  twice 
the  expected  width  extending  in  the  direction  of  the  expected  straight 
line.  (This  description  shows  that,  for  the  purpose,  the  testing  will 
be  done  with  a  rather  narrow  rectangular  test  area,  which  must 
eventually  be  rotated  if  one  tries  to  scan  for  straight  lines  of  differ¬ 
ent  direction. )  Eq  (21)  applies  immediately  and  allows  us  to  determine 
the  limit  where  a  pattern  can  still  be  recognized. 

One  may  ask,  what  will  happen  if  one  tests  for  a  pattern  which 
differs  from  the  pattern  which  occurs  in  the  picture.  Let  us  assume 
that  the  line  where  the  darkness  occurs  has  actually  one  half  the  width 
of  what  we  expect  it  to  have  and  that  its  density  is  twice  what  we 
expect  it  to  have,  i.  e.  that  the  value  a^  is  now  twice  the  value  of  the 
previous  case.  Then  since  the  number  of  grains  in  A  and  in  the 
pattern  is  the  same  as  before,  one  obtains  from  the  calculations  the 
same  value  a^  as  before.  In  reality  a^  is  twice  this  calculated  value. 
The  variance  has  the  same  value  as  before.  Furthermore,  if  the 
proper  width  had  been  used,  then  N'  would  be  one -half  of  the  original 
value  and  the  true  value  of  a^  now  is  \[2  times  the  calculated  value 
of  a^.  This  shows  that  the  scanning  should  be  adapted  to  the  pattern 
as  well  as  possible. 
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A  step  in  p  in  the  x  direction  may  be  described  by 


p  =  -1  for  x  <  0 
p  j  =  +1  for  x  >  0 

To  scan  for  such  a  step  one  must  choose  a  function  6  which  will 
not  register  a  constant  background  density,  or  a  density  gradient, 
or  a  step  in  the  y  direction. 

Let 


P5  =  -1  for  y<0 

p  =  +1  for  y>0 
b 

Then  one  has 

e=  p,  +  X/2+  )'3p3  +  x4p4  +\5P5 

with  the  additional  conditions 

//  0  p  dxdy  ^  0  (22a) 

A  1 


V 


31 


i  =  2,  3  . 


5 


(22b) 


/  J  0  p.  dxdy  =  0  , 


These  are  4  conditions  determining  \  ...  X  .  All  X's  except 

«  D 

X.  turn  out  to  be  zero.  This  is  seen  if  one  assumes  6  to  have  the 
form 

0  =  +  X  p3  (23) 

Since  and  p^  are  odd  functions  of  x  while  p^  ,  p^  and  p,_ 
are  even  functions  of  x,  it  is  immediately  obvious  that  condition 
(22b),  is  fulfilled  by  (23). 

For  an  interval  -a<  x  <  a  ,  one  obtains  from  (22c) 

+  a 

/  0  x  dx  =  0 
-a 

and  hence 


2  a 


In  determining  the  mean  value  we  can  omit  all  contributions  to 
P  which  is  due  to  p^  •  •  •  p^  •  Let 


p  =  p0u  +  »,p1) 

The  average  measurement  is 

+ab  3 

my  =  kaipo 11  pi{pi-  2  i)dxdy 

-a-b 

.  v  ajN 

=  kaip0  ab  *  -j- 
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In  determining  the  variance  it  is  assumed  that  p  is  nearly  constant; 
actually  this  assumption  is  not  needed  with  the  present  choice  of  pg . 
One  obtains. 

a  b  2 

Vv  =  2kpo  /  /  <  1  “  T  J  )  tody 

Y  0-b  2  a 

N 

4 


The  discussion  carried  out  in  conjunction  with  (21)  applies 
again.  In  order  to  obtain  the  same  limit  where  a  pattern  can  still 
be  recognized,  the  number  N  in  the  present  case  must  be  4(1-  y^) 
times  as  large  as  before. 

Quite  obviously,  if  the  step  in  the  photograph  occurs  at  a 
position  different  from  x  =  0  and  we  test  for  a  step  at  x  =  0,  we 
shall  obtain  a  result  which  indicates  the  presence  of  a  step.  The 
same  applies  if  the  step  does  not  extend  parallet  to  the  y-axis,  but 
is  inclined  to  it.  In  principle,  it  might  be  feasible  to  add  conditions 
which  would  guard  against  misinterpretations  of  this  kind,  but 
because  of  the  large  number  of  grains  required,  the  practical  value 
of  such  a  method  is  doubtful.  One  would  probably  make  the  test 
indicated  here  for  different  locations  of  the  jump  and  for  different 
directions  of  it  and  then  place  the  jump  in  the  position  where  the 
indication  is  maximum. 


Recognition  of  Darker  Spots  Against  a  Uniform  Background 


The  previous  section  dealt  with  the  problem  of  measuring  certain 
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defined  quantities  and  of  determining  confidence  limits  for  these 
measurements.  In  this  section  we  treat  a  case  where  we  know 
that  the  light  density  can  assume  two  values,  that  of  the  background 
and  that  of  the  signal  which  raises  the  light  density  slightly  above 
the  background.  We  want  to  discriminate  between  areas  of  this 
background  light  density  and  of  the  higher  light  density. 

To  be  specific,  let  us  assume  that  the  plate  is  divided  into  a 
rather  large  number  of  areas  of  equal  size  (e.  g.  squares)  and 
that  these  squares  have  been  exposed  to  the  background  light  or  to 
the  background  light  plus  the  light  of  the  signal.  In  a  fraction  |3 
of  all  squares  we  have  the  background  light  and  corresponding  to  it 
a  value  as  expected  number  of  grains;  in  the  remaining  fraction 
(1  -  (J)  we  have  background  plus  signal  and  correspondingly  an 
expected  number  of  grains  .  For  convenience  of  discussion,  if 
an  area  has  an  expected  number  of  grains  equal  to  N^,  we  shall  say 
it  "belongs  to  N^".  Let  N^<  .  The  regions  corresponding  to 

the  signal  (N^)  are  to  be  found.  For  this  purpose  we  count  the 
number  of  grains  in  each  area  and  consider  the  areas  below  a  certain 
cut  off  count  m  as  belonging  to  N^,  all  areas  above  the  cut  off  count 
as  belonging  to  N^,  . 

m  should  be  chosen  in  such  a  manner  that  the  probable  number 
of  wrong  decisions  is  minimized.  (This  formulation  is  not  the  only 
possible  one,  one  might  e.  g.  count  an  error  committed  by  missing  a 
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dark  spot  more  heavily  than  an  error  where  a  light  spot  is  considered 
as  dark. ) 

The  probability  of  finding  n  grains  in  an  area  which  belongs  to 
an  expected  value  is  given  by  a  Poisson  distribution.  Therefore, 
if  m  is  the  cutoff  count,  the  probability  of  misjudging  areas  belonging 
to  N  as  areas  belonging  to  is  given  by 


f>  e-N‘  l  Hi 

m  n! 


The  probability  of  misjudging  areas  belonging  to  as  areas  belong¬ 
ing  to  Nj  is 


-N,  m“ i  n2 

(1-P)  e  "2  Z 

n! 


Thus  the  probauili+y  of  making  an  error  is 


Pe=Pe 


Ni  oo  N"  _N2  m-1  Mn 

■  1  S  jl1  +  (1-p)  e  L  S  N2 


ni  n ! 


To  find  a  condition  for  m,  we  increase  and  decease  the  cutoff  count 
m  by  one;  in  either  case  the  error  probability  must  rise.  In  the  first 
case  the  error  probability  changes  by 


r™  _n2  n? 

1  +  (1-P)  e  —4 

n  i  m! 
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In  practice  it  may  be  possible  to  change  the  size  of  the  test  areas 
while  the  ratio  N^/N^  is  fixed.  Of  particular  interest  is  the  case 
where  N^/N^  is  close  to  one.  Let 


Then 


The  expression  x'  does  not  depend  explicitly  upon  or  > 
only  upon  their  ratio.  It  gives  a  constant  shift  in  the  cutoff  count 
which  is  independent  of  the  size  of  the  test  areas  (which  is  charac¬ 
terized  by  Nj  or  ).  The  quantity  x"  is  best  described  in  an 
intuitive  manner  by  the  expression 


6  = 


x"-N 

VN 


i.e.  by  the  quotient  of  x"  -  and 


Nj  .  One  finds 


This  quantity  is  shown  in  Fig.  3;  it  is  rather  close  to  2  • 

Of  prime  interest  are  actual  values  for  the  probability  of  making 
errors.  To  compute  these  the  Poisson  distributions  in  (24)  are  replaced 

♦The  symbol  In  represents  'the  natural  logarithm  of" 
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by  the  corresponding  normal  distributions,  and  the  sums  by  integrals, 
(the  admissibility  of  this  approximation  is  shown  in  Appendix  II).  One 
then  finds 


Pe  =  P 


'2irN,  m 


V 


-(x+|-Ni)^ 

2N"i 


dx+(l-p) 


1 


*mar. 


m 

J 


-(x+i-N,)2 


2N, 


Introducing  in  the  first  integral 


x+i-Nj 

•sTNj 


and  in  the  second  integral 


one  obtains 


P  = 
e 


.  ^  Ui  2 

P  /  e  du  +  (1-p)  /  2e-ii 


u. 


du 


where  u^  =  x1 


1 


1 

+  5  (Yj-IJNj 


u2  =  x 


II  I  1 

N  -  (1-5)  (y-  1)  Yi  2  N. 


The  values  for  these  integrals  in  terms  of  the  upper  limit  can  be 

found  in  tables.  Some  examples  are  given  in  Figs.  4  through  12. 

N  2 

In  the  examples  y .  =  —  =  1.1#  1.5,  and  2.  0.  For  f)the  values 

1  1 

0.  1,  0.  5,  and  0.9  have  been  chosen.  The  value  of  N.  varies  from 
50  to  1000.  The  probability  of  making  an  error  should  be  compared 
with  (1-{J)  which  indicates  in  how  many  of  the  cases  a  signal  is 
present.  As  f}  is  reduced  this  ratio  becomes  rather  high  although 
the  probability  of  making  an  error  is  reduced. 

7.  A  Comparison  of  Two  Techniques  of  Assigning  a  Light  Density 
to  a  Single  Area  of  Lesser  Intensity  than  the  Surroundings 
Assume  now  that  we  look  at  a  given  area,  a  portion  of  the  photo¬ 
graphic  plate,  which  we  shall  designate  as  the  "test  area".  The 
number  of  grains  which  we  count  is  n.  We  want  to  decide  whether 
this  test  area  corresponds  to  an  expected  number  or  to  some 
alternative  expected  number  N^,  N^>N^  (without  specifying  what 
N2  is). 

Two  different  strategies  for  making  this  decision  will  be  com¬ 
pared:  We  choose  a  number  r,  r>l.  In  Strategy  1,  we  count  the 
number  of  grains  in  the  test  area.  If  it  is  below  rN^  we  consider 
this  area  as  belonging  to  N^.  In  Strategy  2,  the  test  area  is  sub¬ 
divided  into  smaller  areas  (for  the  sake  of  argument,  into  nine  sub- 
areas,  let  us  say).  Let  N'^  =  N^/9  .  We  choose  another  number 
r',  r'>l.  Then  we  determine  for  each  of  the  smaller  areas  whether 
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r 


the  count  is  smaller  than  r'  N'^ .  To  the  original  test  area  we  then 
assign  the  count  N,  if  more  than  half  of  the  subareas  have  a  count 
smaller  than  r1  N* . 

Several  questions  immediately  arise: 

(1)  How  shall  we  choose  r  and  r'  ? 

(2)  Which  strategy  is  superior  ? 

(3)  How  much  better  quantitatively  is  the  superior  method  ? 

Intuitively  one  might  consider  Strategy  1  (using  the  whole  area) 

as  superior;  however  Strategy  2  (by  subareas)  may  recommend  itself 
for  practical  reasons.  The  object  of  this  section  is  to  determine 
quantitatively  r,  r'  and  the  probabilities  of  error  in  either  case. 

As  stated  in  Section  6  (cf  p.  35),  there  are  two  possibilities  for 
making  a  wrong  decision: 

(1)  misjudging  the  test  area  as  not  belonging  to  when  it 
actually  does  belong  to  (This  is  called  "Type  I  error"); 

(2)  misjudging  the  test  area  as  belonging  to  when  it  actually 
belongs  to  some  other  number  N^,  different  from 
(This  is  called  "Type  II  error"  ).  In  this  section  is 
greater  than  . 

For  a  number  N^(or  equivalently,  N'^)  and  a  fixed  value  of 
r(orr')  one  can  calculate  ihe  probability  of  making  a  wrong  decision  of 
Type  1  for  each  of  the  strategies.  The  probability  of  making  a  wrong 
decision  of  Type  2  is  then  a  function  of  for  each  of  the  strategies. 


40 


Intuitively  one  can  see  that  the  closer  is  to  the  greater  is 
the  probability  of  making  a  wrong  decision. 

Consider  Strategy  1  (whole  area  method).  We  assume  that  the 
energy  density  of  the  incoming  light  is  constant  over  the  test  area. 

To  fix  the  argument  let  us  temporarily  assume  that  the  test  area  has 
an  expected  value  of  N^.  Under  Strategy  1  we  will  decide  that  the  test 
area  does  not  belong  to  if  the  number  of  grains  counted  is  greater 
than  rN^.  Let  P  be  the  probability  that  we  misjudge  the  test  area 
(with  expected  value  N^)  as  not  belonging  to  N  ,  using  Strategy  1. 

Then 


P  = 


N,  _ 
-  i  2 

e  _  2 


[rNj+1] 


(N^ 

n! 


n 


where  the  brackets  indicate  "largest  integer  in"  . 


Then 


1  r  00 

J  e 


^  ZttNj  rN 


(x-Nj) 

1^7  dx  for  ni  - 5 


Let  u  = 


x-N  j 

n/n[ 


Then  P  is  approximately  given  by 


u/?  i  oo  7 

-  -JL  /  e  "  du  =  -==.  I  -  u  h 
sfT*  ^  2ir  3(r-l)  n/n-J 


du 
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Column  4  of  Table  1  gives  values  for  P  corresponding  to 

(r-l)N/Ny  =  0(.l)  1.2. 

Now  consider  Strategy  2  (subarea  method).  If  the  expected 

number  of  grains  in  the  test  area  is  Nj,  then  the  expected 

number  of  grains  in  each  subarea  is  N'^  where  N'^  =N^/9,  and 

the  probability  of  finding  n  grains  in  any  subarea  is  given  by  a 

Poisson  distribution  with  X.  =  N'^  .  Under  Strategy  2  we  will  decide 

that  a  subarea  does  not  belong  to  N'^  if  the  number  of  grains  counted 

is  greater  than  r'N'  .  Let  p  be  the  probability  that  we  misjudge 

i  s 

a  subarea  (with  expected  value  N'^  )  as  not  belonging  to  N'^  .  The 
subscript  s  denotes  the  subarea  method. 

Then 


N' 


Ps  =  C 


(Njf 

[r'N^+l]  ~n\ 


oo 

E 


N  1  r'N,' 


-  -  (x-N',)2 

e  -  dx  for  N'  >  5 

2N'  1 

1  1 


Let  u  = 


x-N 

n/n! 


Then 


p-  =  ^  ; 


2, 

-  u  f7 

e  * 


(r’-ljN/N^ 


du 


In  the  case  r'  =r,  Column  2  of  Table  I  gives  values  of  p  for 

s 
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(r-1)  */Nt=  0(.l)  1.2. 

Further,  we  will  decide  that  the  original  test  area  does  not 

belong  to  Nj  if  more  than  half  of  the  subareas  have  a  count 

greater  than  r'N'  .  Let  P  be  the  probability  (using  the  subarea 
X  s 

method)  that  we  misjudge  the  test  area  (with  expected  value  ) 
as  not  belonging  to  .  Since  the  number  of  grains  counted  in 
each  subarea  is  independent  of  that  for  any  other  subarea,  this 
probability  is  a  sum  of  binomial  terms,  viz: 


P8  =  Ps  +  (1}  Ps  (1-ps)+  {z)pl  U-Ps^-mIipJ  <UPs) 


9 

where  p  is  the  probability  that  all  nine  subareas  with  expected 
s 

value  N'j  are  misjudged  as  not  belonging  to  N'^  ; 

9  8 

(  i  )  P  (l-p  )  is  the  probability  that  eight  subareas  are  mis- 
1  s  s 

judged  as  not  belonging  to  N'j  (and  one  subarea  is  correctly 
judged  as  belonging  to  N'^); 


(  4  )  Pg  (1_PS)  is  the  probability  that  five  subareas  are  mis¬ 
judged  as  belonging  to  Nj  (and  the  other  four  subareas  are  correctly 
judged  as  belonging  to  N'^). 

There  are  several  possible  ways  to  obtain  P  once  p  is  known. 

s  s 

One  is  to  interpolate  between  values  given  in  a  table  such  as  the 
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Tables  of  the  Cumulative  Binomial  Probability  Distribution  (Harvard 
University,  Cambridge  University  Press,  1955).  Another  possibility 
is  to  interpolate  in  a  table  such  as  Tables  of  the  Incomplete  Beta- 
function  Ratio  (edited  by  Karl  Pearson,  Biometrika  Office,  University 
College,  London  1934)  since  it  can  be  shown  (see  Kendall  [1],  p.  120) 
that 

P  =  I  (5,5)  =  1  -I,  (5,5) 

8  Pfl  UPS 

where  I  (p,q)  is  the  incomplete  Beta-function  ratio.  However,  here 

it  was  more  convenient  to  calculate  P  directly  from  p  with  the  use 

s  s 

of  an  E-101  Burroughs  electronic  computer.  The  results  are  given 

in  Tablel.  For  example,  when  r'=r,  Column  3  of  TabL.  I  gives 

values  of  P  corresponding  to  (r-l)^N' =  0  (0.  1)  1.  2  . 
s 

The  values  of  P  and  p  can  be  obtained  by  using  a  table  such  as 

s 

Tables  of  Probability  Functions,  Volume  II  (Mathematical  Tables 

Projects,  National  Bureau  of  Standards,  1942). 

For  the  case  that  r'  =  r,  the  last  column  of  Table  I  compares 

the  Type  I  error  probability  of  Strategy  1  (whole  area  method)  with 

that  of  Strategy  2  (subarea  method)  in  the  form  of  a  ratio,  P/P  . 

s 

Strategy  1  gives  lower  error  probabilities  for  all  values  listed;  the 
difference  is  more  pronounced  if  the  Type  I  error  probability  is  low 
(50%  error  means  no  judgement  at  all).  Let  r  =  1.  05.  This  would 
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correspond  to  a  case  where  we  want  to  discriminate  between  two 
light  densities  which  differ  by  10%  and  where  we  lay  the  cut  off 
grain  count  into  the  middle.  Let  N1  =  100,  which  corresponds  to 
900  points  in  the  larger  square;  then  (r-1)  's/N'  =  .5  and  the  proba¬ 
bility  of  making  a  Type  I  error  is  0.  1096  with  Strategy  2  versus 
.066=  with  Strategy  1. 

Thus  far  we  have  discussed  only  probabilities  of  wrong  decisions 
of  Type  I  for  the  two  strategies  (see  pag^  41).  Consider  now  the  wrong 
decisions  of  Type  II  for  the  two  strategies,  i.  e.  misjudging  the  test  area 
with  expected  value  N^,  N2  >N^,  as  belonging  to  N^. 

For  Strategy  1  (whole  area  method),  we  now  consider  the  test  area 
to  have  an  expected  value  N^,  where  N2  =  YjNjYj>1.  Under  Strategy  1 
we  will  decide  that  a  subarea  lelor.gs  to  if  the  number  of  grains 
counted  is  less  than  rN^.  Let  Q  be  the  probability  that  we  misjudge  the 
test  area  (with  expected  value  N^)  as  belonging  to  N^. 


Then 


_N2  [rN i]  (N2)n 

Q  =  e  2  - 

0  n! 


1  rN  . 

■ -  J  1  e 

\THn2  -« 


(»-Nz»2 

2N2  dx  for  N2  2  5 


x-N2 
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i  <r"  Y,)^ 

Then  Q  - - /  _  2 

^  -JT,  *'a  n  du 


where  =  y,N ^ 

For  Strategy  2,  a  subarea  has  an  expected  number  of  grains 
where  N'  =y,N  /9»  and  the  probability  of  finding  n  grains  in  any  sub- 

L  XX 

area  is  a  Poisson  distribution  with  X.  =N^.  Under  Strategy  2  we  will 
decide  that  a  subarea  belongs  to  if  the  number  of  grains  counted 
is  less  than  r'N^  .  Let  be  the  probability  that  we  misjudge  a  sub- 
area  with  expected  value  as  belonging  to  N'^. 


Then 


-N.[r'N1]  (N4)” 

%  s0  -sr 


j.r'Ni  (x-Np' 


f2¥Nj  "* 


for  N'>  5 


x-N'z 


Then 


-L-  rul  e~—  du  =  ■  ness-  jU2  -u2/2  j 

TZff.i  2  -«  e  du 

r  (r'-yjN/Nj' 

where  u,  =  - — z —  *  u,  =  - - - 

1  -s/Nj  2  3^ 


r'  N*  -N1 
11  2 


Furthermore,  we  will  decide  that  the  original  test  area  belongs  to 
if  more  than  half  of  the  subareas  have  a  count  less  than  r'N'^.  Let 
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Q  be  the  probability  (using  the  subarea  method)  that  we  misjudge 
s 

the  test  area  (with  expected  value  N^)  as  belonging  to  N^. 

Then,  this  probability  is  again  ~  sum  of  binomial  terms: 


=  <!>  (»-<!,)+<’)  q’  (i-q/+<3»  %  ‘‘-S'3 


9  5  4 

+«I>  is(1-V 


where  q  is  the  probability  that  all  nine  subareas  with  expected 
s 

value  are  misjudged  as  belonging  to  ; 

9  8 

(  )  q  (1-q  )  is  the  probability  that  eight  subareas  with  expected 
Is  s 

value  are  misjudged  as  belonging  to  (and  one  subarea  is 
correctly  judged  as  not  belonging  to  N'^); 


9  5  4 

(  .)  q  (1-q  )  is  the  probability  that  five  subareas  are  misjudged  as 
*  s  s 

belonging  to  NJ  (and  four  subareas  are  correctly  judged  as  not 


belonging  to  N', ). 

X 

One  again  has  the  relationship: 


Q  =  I  (5,5)  =  1  -  I  (5,5) 
8 


where  I^tp.q)  is  the  incomplete  Beta-function  ratio. 

We  now  have  the  necessary  formulas  to  enable  us  to  compare  the 


effectiveness  of  the  two  methods  for  deciding  whether  the  light  density 
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of  the  test  area  is  N^.  To  make  the  discussion  more  explicit  let  us 

fix  the  value  of  y^  at  y^  =  1»  1.  We  set  the  probabilities  of  errors  of 

Type  I,  P  and  P  ,  at  some  common  small  positive  value  (less  than  one). 

8 

The  values  .10,  .05,  and  .01  are  used  here.  Then  we  solve  for  the 

corresponding  values  (r-l)VNj  and  (r'-lJ'sTN^  .  (For  a  given  value  of 

Nj,  we  can  then  find  the  appropriate  r  and  r*.  Values  of  i-  and  r1 

are  given  in  Tables  11(a)  and  11(b),  respectively,  for  P  =P  =.10,  .05,  .01, 

s 

with  values  of  ranging  from  5  to  1000  for  r*. 

With  these  values  of  r  and  r*,  and  a  value  of  y^,  we  then  calculate 

the  probabilities  of  making  Type  II  errors,  Q  and  Q  ,  for  the  two  methods 

s 

(using  the  above  formulas).  Values  of  Q  and  are  listed  in  Table  IQ 
for  P  =  P  =.10,  .05,  .01,  with  values  of  y.  ranging  from  1.1  to  2.0 

8  X 

and  values  of  ranging  from  50  to  1000. 

Note  that  Q  is  always  less  than  Q  for  these  values  of  a.  This 

s 

evidence  tells  us  that  Strategy  1  using  the  whole  area  is  the  better 
method. 

As  an  illustration  of  a  use  of  Tables  Ila  and  b  and  Table  III, 
consider  y^  =  1.  1.  This  corresponds  to  the  case  where  the  back¬ 
ground  light  intensity  is  10%  darker.  If  we  set  r  =  1.  05  (i.  e. ,  we 
lay  the  cutoff  grain  count  for  the  whole -area  method  in  the  middle) 
and  P  =  .  05  then  the  table  for  r  yields  N^  =  1000.  And  the  corre¬ 
sponding  value  of  r1  (the  cutoff  grain  count  for  the  subarea  method) 
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would  be  1.  06  for  P  =  .  05.  Thus,  if  the  probability  of  Type  I  error 

s 

is  set  at  .  05  for  both  methods,  the  probability  of  making  Type  II 


error  is  .  1840  for  Strategy  2  as  compared  with  .0740  for  Strategy  1. 

One  can  proceed  further,  using  the  formulas  for  P  and  Q  (for 
the  whole  area  method),  for  example,  in  the  following  manner: 

Set  P  and  Q  at  pre -assigned  small  values.  Then  one 
obtains  a  relationship  between  y^  and  . 


For  the  case  P  =  Q  =  .05  one  has: 


(r-  1)^  =  1.64485 


and 


r- 


_ L  n/nT  =  -  1.  64485 

*7 


Eliminating  r  from  the  above  two  equations  and  solving  for  y^  one 


obtains 


Yi  = 


<1+ 


nTn 


1 


Figure  13  gives  the  graph  of  vs.  y^  using  the  above  equation. 

as  plotted  corresponds  to  the  minimum  expected  number  of  grains 
for  the  test  area  which  is  required  such  that  neither  of  the  probability 
of  Type  I  and  Type  II  errors  exceeds  .  05  (for  the  whole  area  method). 
For  example,  consider  y^=2.0.  This  corresponds  to  the  situation 


r 


where  the  surroundings  are  twice  as  dark  as  the  test  area.  The 
graph  yields  16  as  the  minimum  required  value  for  the  expected 
number  of  grains  in  the  test  area.  The  values  of  correspond¬ 
ing  to  various  values  of  from  5  to  1000  is  given  in  the  table 
below. 


5  3.  012 

10  2.311 

20  1.  871 

30  1.  691 

40  1.  588 

50  1.519 

100  1.  356 

200  1.246 

500  1.  153 

1000  1.  107 


Another  possible  application  of  the  above  formula  for  P  is  the 
determination  of  a  100(1 -P)%  confidence  lower  limit  for  Nj.  This 
can  be  obtained  by  a  plot  of  rN^  vs  corresponding  to  a  pre¬ 
assigned  value  of  P.  Then  for  a  given  value  of  n,  the  observed 

number  of  grains  (taken  on  the  rN.  axis), one  obtains  L  ,  a  lower 

N1 

100(1-P)%  confidence  limit  for  (taken  on  the  axis.)  Figure 

14  gives  the  graph  of  L  vs  n  for  P=  .  05  (95%  confidence).  To 

1 

illustrate  the  use  of  such  a  graph  suppose  that  10  points  are  counted 
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'I 

in  the  test  area,  i.  e.  n=10.  For  P  =  .05,  the  corresponding  value 

for  is  6.  Hence  a  lower  95%  confidence  limit  for  the  expected 

Ni 

number  of  grains  in  the  test  area  is  6.  In  layman's  language,  this 
means  that  if  we  count  10  points  in  the  test  area  we  can  expect  (with 
a  95%  chance  of  being  right  in  the  long  run)  the  average  number  of 
grains  in  the  test  area  to  be  greater  than  or  equal  to  6. 

If  so  desired,  one  may  obtain  upper  1Q0(1-P)%  confidence  limits 
for  Nj  by  considering  <  1,  and  using  a  technique  analogous  to 
that  described  above  in  this  section. 

8.  Assigning  a  Light  Density  to  a  Single  Area  when  Its  Surroundings 
Contain  Both  Lesser  and  Greater  Intensities 

As  in  Section  7,  let  the  number  of  grains  which  we  count  in  the 
test  area  be  n.  We  wish  to  decide  whther  this  test  area  corresponds  to 
an  expected  number  or  to  an  alternative  number  (without  speci¬ 
fying  N^).  However,  now  the  only  restriction  on  is  that  it  is 
different  from  . 

The  two  strategies  stated  in  Section  7  with  suitable  modifications 
may  again  be  compared.  Choose  two  numbers  r^,  r^  with  rj<l,r  >1. 

In  Strategy  1  (using  the  whole  area)  we  count  the  number  of  grains  in 
the  test  area.  If  it  is  greater  than  r^N^  and  less  than  r^N^  we  con_ 
sider  this  area  as  belonging  to  N^.  Otherwise,  we  decide  that  the 
test  area  does  not  belong  to  N^.  For  Strategy  2,  suppose  the  test  area 
is  subdivided  into  nine  subareas  again.  Choose  two  numbers  r ^ ,  r^ 
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r 


with  r 'i  <  1 ,  r^>l.  Then  we  determine  for  each  of  the  subareas 
whether  the  count  is  between  r'^Nj/9  and  r^  N^/9.  To  the  origi¬ 
nal  test  area  we  then  assign  the  count  if  more  than  half  of  the 
subareas  have  a  count  between  rl  N,/9  and  rl  N./9.  One  now 

11  i.  1 

has  four  numbers  to  choose:  r^,  r^,  r'^,  r^  . 

As  in  Section  7,  it  will  be  seen  that  Strategy  1  (using  the  whole 
area)  is  superior  to  Strategy  2.  We  follow  the  same  development 
as  given  in  Section  7  and  obtain  formulas  for  the  Type  I  errors  P,  P 

s 

and  the  Type  II  errors  Q,Qg. 

Thus  if  P  is  the  probability,  using  Strategy  1,  that  we  misjudge 
the  test  area  (with  expected  value  N^)  as  not  belonging  to  N^,  then 


P  =  e 


-N, 


[riNj 


<Nl>"  + 


00 

2 


V 


(Nt) 

n  =  0  n!  nsfr^Nj+l]  n! 


1 


or  P«-J=-  /rl  Nl.  Ot-Ni)2  ^  + 

v2ttN,  _«o  e  2N ,  's/2irN 


1 


(x-Nj) 

e  ~  dx 


1 


1  r  N 
2  1 


2N 


1 


1  (ViWTr1  H2 

~  1  -  — —  /  e  2  du  for  N  >5 

***  (r1-l)\TN1 

It  can  be  shown  that  good  choices  for  pairs  of  values  of  r^,  r^ 
are  ones  which  yield  symmetrical  limits  in  the  above  integration,  i.  e. 
such  that  (r^-1)  =  -(r^-l).  We  will  consider  only  symmetrical 
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1 


limits  here 


(r  -l)^N, 


Thus 


2  '21 
1  '  tBT  !0  e'f  du 


Table  IV  gives  values  of  r^  and  (for  various  values  of  Nj)  such 
that  P  =  .10,  .05,  .01. 

Now  if  Q  is  the  probability,  using  Strategy  1,  that  we  misjudge 
the  test  area  (with  expected  value  N^)  as  belonging  to  N^,  then 


Q  =  e 


N  [^N^ 
i"  1  2 

[ r  jNj+1] 


(N2)n 


or  Q 


n/2^N 


z  ,lNl 


f'?.Ni  J*-N2) 

J  e  2N?  dx  = 


_j_  fr?N1  Jx-N2)  1  Jr2-ylWNl 

Q~^f?2  \  N6  *=  3R  f  e 

1  1  J  - 

(r1-Y1)'/N1 

\/y, 

3{C  1 

v/here  N^  =  YjNj  ,  y  ^  >  0,  and  >  5 


Substituting  for  r  ^  in  terms  of  r^  [using  the  symmetrical  relation 
(r^-1)  =  - (r^-l)  ]  we  obtain 


(r2-l)N/N'1  _  (yj-I)^ 

Q  w  L»  /  "sTvT  \[y\  -  u2  - 

n/2u  1  _  i  e  -y  du  for  y^N  >5 

-(r2-l)'s/N1  _  (yj-D-y/Nj 

^  ^yT 

♦Note  that  y^  may  be  less  than  or  greater  than  one,  since  N2  may  be 
greater  than  or  less  than  Nj  for  this  case. 


53 


r 


As  before,  Pg  is  the  probability,  using  Strategy  2,  that  we 
misjudge  the  test  area  (with  expected  value  N^)  as  not  belonging 
to  Nj.  By  analogy  to  the  argument  used  in  the  previous  section, 
we  obtain  the  following  formula  for  this  case: 

Pc  =  Xn  (5'5) 

S  Ps 

where  I  (p,  q)  is  the  incomplete  Beta-function  ratio 

(r-  - 

2  ^  2  1 

and  p  ~  1  -  —r—  /  e"H_  du  for  N‘ >  5 
s  v2ir  0  ^  1 

As  before,  p  is  the  probability  that  we  misjudge  a  subarea  (with 

expected  value  N'^)  as  not  belonging  to  N'^  . 

Finally,  if  Q  is  the  probability,  using  Strategy  2,  that  we  mis- 
s 

judge  the  test  area  (with  expected  value  N^)  as  belonging  to  N^,  then 

Q  =1  (5,5) 

8  qs 

where  I  (5,5)  is  the  incomplete  Beta  function  ratio 

2 
u 

2  du 


with  =  Yj  Nj  and  ^  5 


and  q  __ 


VTir  1 


(r2-V^'l  J 

^1  e  2  du  = 

(r*  -  Yi)Vn> 

n/y7 


(r2-Vi> 


x  ,  r-  N^N. 

n/2tt  e 

-(rz  +  y/i]  ^ 


1 
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As  before,  q  is  the  probability  that  we  misjudge  a  subarea  (with 
s 

expected  value  Nl,  )  as  belonging  to  N'^  . 

Proceeding  as  we  did  in  Section  7,  we  set  P  and  P  (the  Type  I 

8 

errors  for  the  two  methods)  at  some  common  value.  Again  the  values 

.10,  .05,  and  .01  are  used  here.  Then  we  obtain  the  corresponding 

values  of  (r  -1)  and  ( r '  - 1 )  V  N using  the  formulas  for  P  and  P 

£  1  £  1  8 

above.  For  a  given  value  of  N^,  we  can  then  find  the  appropriate  r^ 
and  r^  .  Furthermore,  we  cap  obtain  the  appropriate  r^  and  r^  using 
the  symmetry  relations. 

Values  of  r^,  r^  are  given  in  Table  IV(a)  and  values  of  r^,  r^  are 

given  in  Table  IV(b)  for  P=P  =.10,  .05,  .01  with  values  of  N  ranging 

s  1 

from  5  to  1000  for  (r  ,  r  )  and  from  50  to  1000  for  ( r' ,  r'). 

1  L  1  b 

For  given  values  of  r  and  r',  and  a  value  of  ,  we  can  calculate 

Q  and  Q  (the  probabilities  of  Type  II  errors  for  the  two  methods), 
s 

These  are  given  in  Table  V  for  P=P  =.10,  .05,  .01  with  values  of  y  . 

8  * 

ranging  from  0..5  to  2.  0  and  values  of  ranging  from  50  to  1000. 

As  in  the  previous  section,  an  examination  of  Table  V  reveals  that 

Q  is  always  less  than  Q  for  the  values  considered.  This  evidence 

s 

tells  us  that  Strategy  1  (using  the  whole  area)  is,  as  expected,  the  better 
method.  For  example,  consider  y^=0.  9.  This  corresponds  to  the  case 
where  the  test  area  is  10%  lighter  than  the  background  light  intensity. 

If  we  set  r^  =  .94,  r^  =  1.  06  and  P  =  .05  then  N ^  =  1000.  The  corre¬ 
sponding  values  of  r'j,  r^  (for  Pg  =  .05)are  r^  =  .89  and  rj,  =  l.ll.  Thus, 
if  the  probability  of  Type  I  error  is  set  at  .05  for  both  methods,  the 
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probability  of  making  Type  II  error  is  .10  for  the  whole  area  method 
as  contrasted  with  .57  for  the  subarea  method. 

Again,  one  can  set  P  and  Q  at  preassigned  small  values  and 
obtain  a  relationship  between  y  ^  and  N^.  For  the  case  P  =  Q  =  .05 
this  expression  is  different  from  that  obtained  in  Section  7. 

Thus,  when  P  =  Q  =  .05  one  has: 


(r2  -  1  )<sTN  =  1.95996 


r,-y.  + 

Also  _£ — LnTNj  1.64485  for  Yj  >  1 


("r2  "  Yi  +  2)nTn  + 

and  - 7 -  »  1.64485  for  y  <  1 

nTYi  1 


Eliminating  r0  from  these  equations  and  solving  for  y  we  obtain 
Yj  =2n/N1[c1  sgn(y1-l)  +  N/N'1]  +  c2  +  c2  sgn^-l)  ^WnJcj  sgri^-lH-s/lT^+c 


2N 


1 


where  c  ^  =  1.  95996 

+ 

c2=  1.  64485 


for  Yj>0 


Figure  15  gives  the  graph  of  vs  y^  corresponding  to  the  above  equa¬ 
tion.  As  in  the  previous  section  this  value  of  Nj  (for  a  given  value  of  y^) 
is  the  minimum  expected  number  of  grains  for  the  test  area  which  is 
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required  such  that  neither  of  the  probabilities  of  Type  I  and  Type  II 
errors  exceeds  .05  (for  the  whole  area  method).  The  values  of 
corresponding  to  various  values  of  from  5  to  1000  is  given  in  the 
table  below. 


N1 

Yl 

(Y1>1) 

Yl 

(Yi<D 

5 

3.  190 

.0198* 

10 

2.  431 

.  1674* 
+ 

20 

1.  952 

.  3455 

30 

1.  756 

.  4424 

40 

1.  643 

.  5052 

50 

1.  569 

.  5503 

100 

1.  390 

.  6694 

200 

1.  270 

.  7600 

500 

1.  167 

.  8447 

1000 

1.  117 

.  8890 

As  an  example,  the  graph  in  Figure  15  yields  =  18  for  y^  =  .31  and 
2.  1  .  Corresponding  values  of  for  these  values  of  are  respec¬ 
tively  5.  6  and  37.  8.  Thus  if  it  is  desired  to  distinguish  between  the 
test  area  and  surrounding  areas  containing  either  0.  31  times  as  many 
grains  per  area  (on  the  average)  as  the  test  area,  or  2.  1  times  as  many 
grains  per  area  (on  the  average)  as  the  test  area,  then  the  minimum 
required  number  of  grains  in  the  test  area  (on  the  average)  is  18  if 
neither  of  the  probabilities,  P  and  Q,  of  Type  I  and  Type  II  errors  is  to 
exceed  .05.  As  another  example  of  the  use  of  Figure  5,  suppose  it  is 

♦These  figures  are  not  accurate  to  four  places  since  <  5  for  these  cases. 
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desired  to  distinguish  between  the  test  area  and  surrounding  areas 
containing  either  0.  5  times  as  many  grains  per  area  or  2.  0  times 
as  many  grains  per  area  as  the  test  area.  Then  =  0.  5  yields 
=  40  and  =2.  0  yields  =  18  in  the  graph.  The  larger  of  the 

two  values  for  N^,  namely  40,  then  is  the  minimum  required  number 
of  grains  in  the  test  area  if  neither  of  the  probabilities,  P  and  Q,  of 
Type  I  and  Type  II  errors  is  to  exceed  .05. 

To  obtain  100 (1-P)%  confidence  intervals  for  N^,  plot  r^N^, 
r^Nj  versus  for  a  pre-assigned  value  of  P.  Then  for  a  given 
value  of  n  (on  the  r  N  ,  r  N  axis)  one  obtains  1 00 (1-P)%  lower  and 

XX  b  1 

upper  confidence  limits  for  N^,  and  (on  the  axis),  from  this 

graph.  Figure  16  gives  the  graphs  of  n  versus  (LN  ,  UN)  for  P=.10 
.05  .01.  For  example,  suppose  that  10  points  are  counted  in  the  test 
area.  Thus  n=l'0.  For  P  =  .05,  the  corresponding  values  for  and 

are  7  and  19  respectively.  This  means  that  if  we  count  10  points 
in  the  test  area  we  can  expect  (with  a  95%  chance  of  being  right  in  the 
long  run)  the  average  number  of  grains  in  the  test  area  to  be  between 
7  and  19. 

Conclusion 

This  report  analyzes  procedures  for  the  evaluation  of  strongly 
enlarged  photographs  by  electronic  scanning  techniques.  The  enlarge¬ 
ment  is  so  great  that  the  grain  --iructure  of  the  photograph  limits  the 
accuracy  of  the  evaluation.  Two  sources  of  uncertainty  are  present 
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in  this  process:  (1)  The  quantity  to  be  determined  for  a  certain  point 
of  the  photograph  (in  many  cases  the  local  energy  density  of  the  elec¬ 
tromagnetic  waves)  must  be  approximated  by  an  integral  over  an  area 
containing  this  point;  (2)  In  the  determination  of  this  integral  an  uncer¬ 
tainty  occurs  because  of  the  randomness  in  the  grain  distribution.  To 
reduce  the  first  uncertainty  the  test  area  should  be  taken  small;  to 
reduce  the  second  one  the  test  area  must  be  taken  large. 

Regarding  the  second  source  of  errors  it  is  shown  how  the  variance 
of  the  results  can  be  expressed  by  a  integral  containing  the  light  density 
and  a  weight  function.  By  a  suitable  choice  of  the  weight  function  the 
variance  is  minimized.  This  procedure  is  illustrated  in  several 
examples.  Included  is  the  problem  of  scanning  for  special  patterns  in 
the  light  density.  The  uncertainty  mentioned  under  (1)  is  connected 
with  the  size  of  the  test  area  and  the  properties  of  the  photograph.  In 
principle  the  combined  variance  due  to  the  two  sources  of  error  can  be 
minimized. 

Further  investigations  are  concerned  with  the  discrimination  between 
two  known  light  densities  and  a  special  stepwise  procedure  in  evaluating 
the  light  density  of  the  test  area. 

Formulae  for  the  confidence  limits  on  the  true  grain  count  are  given 
for  different  cases.  The  confidence  limits  depend  naturally  upon  the 
average  number  of  grains  occurring  in  the  test  area.  (If  several  photo¬ 
graphs  of  the  same  object  exist,  the  number  of  points  in  a  given  area  are 
added).  For  moderate  accuracy  requirements  (5%  to  10%  uncertainty)  a 
lower  limit  is  given  by  1000  grains. 
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APPENDIX  I 

DERIVATION  OF  FORMULAS  FOR  My  AND  Vy 

Uniform  Distribution  of  Grains  in  an  Area 

The  statistical  theory  given  in  this  appendix  is  not  self-contained. 
F*>r  the  convenience  of  the  reader,  references  have  been  made  to  one 
book  only,  ,rThe  Advanced  Theory  of  Statistics,  Vol.  1"  by  M.  G. 
Kendall,  [1], 

Let  us  assume  that  a  constant  number  of  grains,  N,  are  distri¬ 
buted  with  uniform  probability  in  an  area  y  •  Further,  assume  that 
grains  arising  in  some  subarea  Au  of  y  are  statistically  independent 
phenomena.  One  may  proceed  by  raising  the  following  simple  question: 
If  one  counts  u  grains  in  an  area  Aw  what  are  the  limits  of  error  on 
this  number  u? 

Since  the  production  of  grains  occurs  with  uniform  probability, 
it  follows  that 

Aw 

P=—  (Al) 

is  the  probability  that  a  grain  is  found  in  Aw,  provided  that  we  know 
that  it  is  in  y. 

Let  X  represent  the  number  of  grains  counted  in  Aw.  It 
will  vary  according  to  the  particular  photograph  .*  Since  p  is  a 
constant  for  an  area  of  size  Aw  in  y  it  follows  that  the  probability 
that  X  is  equal  to  a  number  u  (written  as  P^(u)),  is  given  by  a 

♦Statisticians  refer  to  it  as  a  "random  variable."  ([1]  page  173) 
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binomial  distribution 


P  (u)=  CN  pU(l-p)N_U  u  =  0,  1 . . .  N  (A2) 

A  11 

For  a  discussion  of  the  binomial  distribution  see  [1]  page 
116.  In  general,  for  our  problem  N  will  be  very  large  and 
p«  1 .  Let  X  =  N  represent  the  average  number  of  grains 
expected  in  Aw.  If  N-»»  and  Aw-*0  in  such  a  manner  that  X 
is  a  constant  less  than  5,  P  (u)  approaches  a  Poisson  distribu- 

A 

tion  with  mean  X  ;  i.  e. 

,  ,  a'  Su 

Px(“)  =  — —  (A3) 

It  is  well  known  that  for  the  Poisson  distribution  the  mean 
and  the  variance  are  both  equal  to  X  .  For  a  discussion  of  this 
distribution  see  [l]  p.  120. 

At  this  point,  it  is  worthwhile  to  introduce  the  notion  of  the 

"characteristic  function"  <}>  (t)  which  is  uniquely  associated  with 

A 

P  (u) .  The  introduction  of  «J>  ( t)  facilitates  greatly  the  deriva- 
A  A 

tion  of  formulas  for  and  (see  page  4  of  this  report  for 

definitions  of  M  and  V  ).  4>  (t)  is  defined  by 

I  i  A 

4>v(t)  =  S  eitU  P  (u)  (A 4) 

x  u  =  0  X 
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<(>  (t)  is  known  in  statistics  as  the  "characteristic  function  of  Pv(u) 
X  X 

{corresponding  to  the  random  variable  X).  See  [l]  page  90  for 
further  discussion  on  this  subject.  Using  (A3)  and  (A4)  we  obtain 


i  .  v  »  00  t  \  ^  V  ^ 

.  .  ltu  -x  u  _  -X  ”  Xe  ) 

<j>v(t)  =  S  e  e _ X  =  e  2  - —  = 

X  u  =  0  u  =  0  u! 


-X  X  e 
e  e 


Me^-l) 


Now  consider  y  as  being  subdivided  into  a  finely  divided 
grid  of  squares  where  each  square  has  the  same  area  Au.  Label 
the  squares  of  the  grid  as  Auj^  .  .  Aw.,  .  Aw^  (Fig.  17). 

We  now  consider  a  random  variable  Y 


Y  = 


n 

2  0.  X. 
i  =  1  1  1 


\A6) 


where  X^  =  number  of  grains  counted  in  Au^ 

0.  =  weight  assigned  to  the  points  in  Au^  . 

It  is  assumed  that  X^...  X^  are  statistically  independent 
of  one  another.  The  corresponding  characteristic  function  of  Y 
is  obtained  as  follows: 

First  one  finds 
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The  X.'s  are  mutually  independent.  Hence  for  a  fixed  set  of 


values  u.,  . . .  u  of  X.  ...  X  respectively  we  have: 
lnln 


u 


Pv(2  0.u.)  =  P(u.,u  ,  ...  ,  u  ) 
»  1  1  l  c  n 


(A  8) 


n 


=  ii  P(u4) 


i=l 


See  [1]  page  21. 
Hence 


.  .  *  «  it(0 ,u  +. . .  +0  u  ) 

<|>  (t)  =  Z  .  .  . .  2  e  11  n  n  n 

Y  u  =0  u  =0  Pv(  Z  0.  u.  ) 

1  n  Y  i  -  1  1  1 


_  ®  “  it  ( 0  u  .  .  .  +  0  u  )  n 

-Z  ,...Z  e  11  nn.irP(u.) 

u  =0  u  =0  1 

1  n 


(A9) 


2  it0 .  u 


=  Z_0  e“~riP(u  )  Z  elt02U2  P(u  )...  Z  ellUnP(u) 

u!  0  1  u2=0  2  u  _0  n; 


itUn 


n 


MoVx.«*> 

li 


see [l]p242 


(A  7)  and  (A9)  yield 


n 


.  /  it0.  “  .  it0.  .. 


(A  10) 
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4»Y^)  =  exP  {“T-.  Au>.  (  eit0i-l )  }-*exp  {  ~  //[  eit0  ^X’  ^ -l  ]dxdy}  (A  14) 
•  i  1  ^  « 


assuming  that  the  double  integral  exists. 
From  this  one  obtains 


log<v<.)=  S//[ei*e<-‘,y>.l]dxdy 


This  yields 


V  y  r=l 


=  2  Mf[2  //  0r(x,  y)  dxdy] 

r  =  l  r!  Y  y 


(A  15) 


Thus 


K  =  £  ;;er(x,y)d*dy 
y  Y 


(A  16) 


M  =  ^  //  0  (  x,  y)  dxdy 
y  Y  Y 


(A  17) 


Vy=  *  //02(x,y)dxdy 


(A  18) 
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General  Case 


It  has  been  assumed  in  the  text  that  the  grains  arise  by  a  random 
process  in  such  a  manner  that  the  probability  that  a  grain  is  produced 
in  a  given  area  is  proportional  to  the  light  flux  across  that  area.  On 
the  basis  of  this  assumption,  it  is  natural  to  make  the  following  defi¬ 
nition  for  p,  where  p  is  the  probability  of  a  grain  in  Aw^  given  that 
the  grain  is  in  y.  Define  p  as: 


//  p (x, y) dxdy 
Ab>i 

p=— - 

J7  p(x,y) dxdy 

Y 

where  p(x,  y)  is  the  density  of  the  light  flux. 

Then,  since  N=k//  p(x,y)dxdy  by  Equation  1 
Y 

p=  ^  //  P(x.y)  dxdy 
Au> 


(A19) 


(A20) 


Note  that  when  p(x,y)  is  a  constant  and  Au.  has  area  A w, 

p  =  -  an(i  this  reduces  to  the  case  discussed  in  Section  1. 

Y 

We  shall  assume  that  the  region  y  has  been  subdivided  by  some 

suitable  method  into  subdivisions  Au.  such  that  p,  as  defined  in 

(A  19),  is  the  same  for  any  Aw^  thus  produced.  Then  the  discussion 

in  Section  1  can  be  applied  again.  Here,  X  =Np=  k  //  p  (x,  y)dxdy  is 

Au>j 

the  average  number  of  grains  expected  in  Aw..  By  definition  this  is 
a  constant  for  any  i. 
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One  obtains  in  a  manner  similar  to  that  given  in  Section  1, 


log  <|>Y  (t)=‘k  £  [  //  p(x,y)dxdy][  elt6l-l] 

i  =  1  Aw. 


K  =  k  E  [  //  p  (x,  y)  dxdy  ]  0* 

1  *  =  1  Aw{ 

If  one  considers  6.  as  a  discontinuous  function  of  x  and 
i 

assumes  in  each  region  Aw.  the  value  0^  defined  above, 
v/rite  the  last  expression 

Kr  =k  J7p(x,y)  0*  (x,  y)  dxdy 

Y 

If  Aw^O,  0^(x,  y)  tends  to  a  function  0  (x,  y),  i.  e. 

[  O.  (x»  y)  -  0  (x,  y)  ]  =  0 

Then 

K  =  k  //  p  {x ,  y)  £ (x,  y)  dxdy 

Y 

Hence 

My  =  ^  //  P  (x»  y)  9  (x,  y)  dxdy 

VY  =  k  //  p  (x,  y)  02(x,  y)  dxdy 

Y 

REFERENCE: 


(A21) 

y  which 
one  can 


(A22) 

(A23) 


1.  Kendall,  M.  G. ,  The  Advanced  Theory  of  Statistics, 
Vol.  I,  1948,  London,  Charles  Griffin  and  Company,  Limited. 


APPENDIX  II 


COMPARISON  OF  POISSON  DISTRIBUTION  AND  NORMAL  DISTRIBUTION 

The  assumption  that  the  grains  of  the  photograph  arise  in  a 
random  fashion  and  that  the  creation  of  a  new  grain  can  be  considered 
as  an  event  independent  of  the  creation  of  all  previous  grains  leads  to 
a  Poisson  distribution.  To  be  specific  if  the  average  number  of  grains 
to  be  found  in  a  certain  test  area  is  N  then  the  probability  of  finding 
in  such  a  test  area  n  grains  is  given  by 

-N  N11 

p(n)  =  e  —  (A24) 


For  large  values  of  n  this  expression  can  be  approximated  by  a 
normal  distribution.  It  is  obtained  if  one  expresses  n!  by  means 
of  Stirling's  formula  and  then  carries  out  certain  simplifications, 
using  the  assumption  that  the  number  "n"  differs  from  the  average 
value  N  by  only  a  small  amount.  One  thus  obtains  p^,  the  proba¬ 
bility  density  of  finding  n  grains  in  a  given  area  of  size  dn(  if  the 
expected  number  is  N)  as  the  normal  law: 


P(n)  =  ‘TOTS 


-  (n-N+ 
2N 


(A25) 


In  Figs.  18-23  the  expressions  (A24)  and  (A25)  are  plotted  for 
N  =  10  (10)  60.  One  recognizes  that  the  normal  distribution  gives 
a  good  approximation  for  N  as  low  as  10  (In  fact  it  is  fairly  good  for 
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N  as  low  as  5).  Thus  it  appears  to  be  justified  if  one  takes  the 
normal  distribution  as  approximation  for  the  Poisson  distribution 
whenever  such  an  approximation  is  needed. 
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Table  I 


Values  of  P  and  P  for  the  case  when  r'=  r  (  y  >  1) 

S  X 


DnTn* 

p 

P 

P 

P/P 

l 

s 

s 

s 

0 

.  5000 

.  5000 

.  5000 

1.  000 

0.  1 

.  4602 

.4028 

.  3821 

.949 

0.2 

•  4207 

.  3113 

.2743 

.881 

0.3 

.  3821 

.2301 

.  1841 

.  800 

0.4 

.  3446 

.  1628 

.  1151 

.  707 

0.5 

.  3085 

.  1097 

.  0668 

.609 

0.6 

.  2743 

.0704 

.  0359 

.  510 

0.7 

.  2420 

.0429 

.  0179 

.416 

0.8 

.  2119 

.0250 

.0082 

.  328 

0.9 

.  1841 

.0138 

.  0035 

.251 

1.0 

.  1587 

.0072 

.0014 

.  187 

1.  1 

.  1357 

.0036 

.0005 

.  134 

1.2 

.  1151 

.0017 

.0002 

.094 
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Table  11(a) 


Values  of  r  for  Yj  >  ^ 


N1 

DU 

P  =  P  =.05 
s 

| 

1 

5 

1.  573 

1.736 

2.040 

10 

1.405 

1.  520 

1.  736 

20 

1.  287 

1.  368 

1.  520 

30 

1.  234 

1.  300 

1.425 

40 

1.  203 

1.260 

1.368 

50 

1.  181 

1.233 

1.  329 

100 

1.  128 

1.  164 

1.233 

200 

1.  091 

1.  116 

1.  164 

500 

1.  057 

1.  074 

1.  104 

1000 

1.  041 

1.  052 

1.  074 

Table  11(b) 

Values  of  r'  for  y^>l 

N. 

p=p  =  .10 

P  =  P  =.  05 

P=P  =.  01 

1 

s 

s 

s 

50 

1.221 

1.284 

1.403 

100 

1.  156 

1.  201 

1.285 

200 

1.  Ill 

1.  142 

1.202 

500 

1.  070 

1.090 

1.  128 

1000 

1.  049 

1.  064 

1.090 

Table  HI 


Values  of  Q  and  Q  for  v.>l 
8  '  1 

P  =  P  =.10 
8 


V 

50 

V 

100 

Nl  = 

200 

Nl  = 

500 

N1 

=  1000 

Y,  Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

*1  8 

8 

8 

8 

8 

1.  1  .7488 

.  7081 

.6706 

.  6058 

.  5469 

.4497 

.2998 

.  1814 

.  1061 

.0365 

1.2  .5449 

.4518 

.  3723 

.2560 

.  1722 

.0790 

.  0151 

.0018 

.  0002 

.0000 

1.  3  .  3445 

.2307 

.  1512 

.  0659 

.  0275 

.  0047 

.0001 

.0000 

.  0000 

.0000 

1.  5  .0938 

.  0329 

.  0110 

.  0012 

.0001 

.  0000 

.  0000 

.0000 

.  0000 

.0000 

2.0  .0000 

.  0000 

.0000 

.0000 

.  0000 

.0000 

.  0000 

.0000 

.  0000 

.0000 

P  =  P  =.05 
8 


Nl  = 

50 

Nl  = 

100 

Nl  = 

200 

Nl  = 

500 

V 

1000 

Y1 

Q 

8 

Q 

Q 

8 

Q 

Q 

s 

Q 

Q 

8 

Q 

Q 

8 

Q 

1.1 

.8458 

.  8144 

.  7852 

.  7307 

.6795 

.  5870 

.4300 

.2865 

.  1840 

.0740 

1.2 

.  6724 

.  5834 

.  5032 

.  3729 

.2702 

.  1400 

.0331 

.  0049 

.  0007 

.  0000 

1.3 

.4682 

.  3380 

.2384 

.  1173 

.0546 

.  0114 

.  0004 

.  0000 

.  0000 

.  0000 

1.5 

.  1538 

.  0613 

.  0231 

.  0031 

.  0004 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

2.0 

.0018 

.  0001 

.0000 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

P  =  P  =.01 
s 


N:  =  50 

II 

zT 

100 

V 

200 

Nl= 

500 

II 

H 

z 

1000 

*1 

Q  Q 

s 

Q 

8 

Q 

Q 

s 

Q 

Q 

8 

Q 

Q 

8 

Q 

1. 1 

.9528  .  9387 

.9257 

.  8970 

.  8692 

.  8078 

.  6847 

.  5343 

.4039 

.2127 

1.2 

.  8588  .  7975 

.  7380 

.  6171 

.  5068 

.  3234 

.  1127 

.  0251 

.  0049 

.0001 

1.3 

.  7002  .  5714 

.4573 

.  2773 

.  1586 

.  0464 

.0029 

.  0001 

.  0000 

.  0000 

1.5 

.  3234  .  1617 

.  0755 

.  0145 

.  0025 

.  0001 

.0000 

.  0000 

.0000 

2.  0 

.0075  .0004 

.0000 

.  0000 

.  0000 

.  0000 

.0000 

.  0000 

.  0000 
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Table  IV(a) 


Values  of 

r  j  and  i 

'2  lor  V1 

P  =  P 
a 

=  .10 

! 

P  =  P 
s 

=  .05 

P  =  P 
s 

=  .01 

N1 

rl 

r2 

rl 

r2 

rl 

r2 

5 

.264 

1.  736 

.  123 

1.  877 

-.  152 

2.  152 

10 

.480 

1.  520 

.  380 

1.620 

.  185 

1.815 

20 

.632 

1.  368 

.  562 

1.438 

.424 

1.  576 

30 

.700 

1.  300 

.  642 

1.  358 

.  530 

1.470 

40 

.740 

1.260 

.690 

1.  310 

.  593 

1.407 

50 

.767 

1.233 

.  723 

1.277 

.636 

1.  364 

100 

.836 

1.  164 

.  804 

1.  196 

.  742 

1.  258 

200 

.  884 

1.  116 

.  861 

1.  139 

.  818 

1.  182 

500 

.926 

1.  074 

.912 

1.088 

.  885 

1.  115 

1000 

.948 

1.052 

.938 

1.062 

.919 

1.  081 

Table  IV(b) 

Values  of 

and  for 

V1 

P  =  P 

8 

P  = 

P  =.05 
s 

P  = 

P  =.01 
s 

Ni 

ri 

r2 

r1 

i 

rz 

r* 

i 

r2 

50 

.561 

1.439 

.  513 

1.487 

.419 

1.  581 

100 

.690 

1.  310 

.  656 

1.  344 

.589 

1.411 

200 

.781 

1.219 

.  757 

1.243 

.710 

1.290 

500 

.  861 

1.  139 

.  846 

1.  154 

.816 

1.  184 

1000 

.902 

1.098 

.891 

1.  109 

.  870 

1.  130 
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Table  V 


Values  of  Q  and  Q  for  v,  ^1 
s  '1 

P  =  P  =  .  10 
s 


N1 

=  50 

Nl  = 

100 

N1 

=  200 

500 

Nx  =  1000 

V. 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

*  1 

s 

s 

s 

s 

s 

0.5 

.  3060 

.  0037 

.  0142 

.  0000 

.  0000 

.  0000 

0.  6 

.  7474 

.  19  it 

.  3982 

.  0193 

.  0526 

.  0001 

.  0000 

.  0000 

0.  7 

.  7992 

.2845 

.  5218 

.0526 

.  1297 

.  0010 

.  0002 

.  0000 

.  0000 

.0000 

0.8 

.  8876 

.  6014 

.7945 

.  3456 

.  5667 

.  0929 

.  1020 

.  0008 

.  0017 

.0000 

0.9 

.9127 

.  8320 

.  8954 

.7490 

.  8561 

.  5954 

.  7108 

.  2666 

.4456 

.0549 

1.  1 

.  8508 

.8019 

.  8316 

.  7248 

.  7904 

.  5853 

.  6534 

.  2864 

.4247 

.0740 

1.2 

.  7574 

.5808 

.  6706 

.  3725 

.4945 

.  1399 

.  1366 

.  0049 

.0078 

.0000 

1.  3 

.  6205 

.  3375 

.4383 

.  1173 

.  1779 

.  0114 

.  0047 

.  0000 

.  0000 

.0000 

1.  5 

.  3047 

.  0613 

.  0860 

.0031 

.  0038 

.  0000 

.  0000 

.  0000 

2.  0 

.0089 

.  0001 

.  0000 

.  0000 

.  0000 

.  0000 

P  =  P  = 
s 

.05 

N1  = 

50 

Nl  = 

100 

Nl  = 

200 

Nl  = 

500 

Nj  =1000 

•y, 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

1 

s 

s 

s 

s 

s 

0.  5 

.4552 

.  0129 

.  0357 

.  0000 

.  0000 

.  0000 

0.6 

.8479 

.  3134 

.  5361 

.  0463 

.  1003 

.  0004 

.  0000 

+ 

+ 

0.7 

.  8847 

.4235 

.  6555 

.  1069 

.2140 

.  0032 

.  0008 

.  0000 

.0000 

0.8 

.9427 

.  7290 

.  8791 

.  4821 

.  6916 

.  1658 

.  1697 

.  0025 

.  0043 

.  0000 

0.9 

.9579 

.9042 

.9472 

.  8433 

.9215 

.  7173 

.  8127 

.  3855 

.  5707 

.  1 02$ 

1.  1 

.9173 

.  8784 

.  9037 

.  8176 

.  8732 

.  6979 

.  7611 

.  3961 

.  5427 

.  1258 

1.2 

.  8471 

.  6898 

.  7749 

.  4853 

.  6113 

.  2139 

.2051 

.  0109 

.  0153 

.  0000 

1.  3 

.  7297 

.4436 

.  5527 

.  1808 

.  2563 

.  0226 

.  0093 

.  0000 

.  0000 

1.  5 

.4047 

.0991 

.  1322 

.  006t 

.  0074 

.  0000 

.  0000 

2.  0 

.0156 

.  0002 

.  0001 

.  0000 

.  0000 
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Table  V  (continued) 


N1 

=  50 

Nl  = 

100 

V 

200 

N1 

=  500 

Nl  = 

1000 

Y1 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

Q 

s 

s 

s 

8 

8 

0.  5 

.7460 

.0874 

.  1504 

.0003 

.  0003 

.0000 

.0000 

0.6 

.9573 

.  6056 

.7792 

.  1768 

.2708 

.0044 

.0004 

.0000 

0.7 

.9703 

.  7065 

.  8574 

.  3061 

.4462 

.  0232 

.  0056 

.0000 

.0000 

0.8 

.9882 

.9030 

.9665 

.  7401 

.  8724 

.  3888 

.  3674 

.  0170 

.0214 

.0000 

0.9 

.9922 

.9753 

.9893 

.9516 

.9812 

.8896 

.9354 

.6399 

.  7829 

.2682 

1.  1 

.9792 

.9617 

.9740 

.  9332 

.9614 

.  8659 

.9037 

.  6270 

.7491 

.2880 

1.2 

.9489 

.  8554 

.9104 

.7004 

.  8002 

.4088 

.  3816 

.  0417 

.  0483 

.0003 

1.3 

.  8817 

.6549 

.7501 

.  3549 

.4433 

.0719 

.  0307 

.  0001 

.  0001 

.0000 

1.5 

.  6050 

.2166 

.2612 

.  0239 

.  0234 

.  0001 

.  0000 

.0000 

.  0000 

2.0 

.0403 

.0007 

.0004 

.0000 

.  0000 

.0000 
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APPENDIX  HI 


I.  Computational  procedure  for  finding  r,  given  a  an  assigned 
value  of  P{x). 

2 

1  .»  _ _ 

P(x,=^?r'*e  2  du 

where  x  =  (r-  Dn/nJ 
Now  set  P(x)  equal  to  a,  i.  e. 

2 

.  _u _ 

P(x)=  — pass-  / *°e”  2  du  =  a  (a) 

Given  a,  the  problem  is:  compute  P  (x).  To  solve  this  problem  on 
Burroughs  E101  Electrodata  Electronic  Computer,  we  proceed  in  the 
following  manner.  We  replace  P(x)  in  (a)  by  the  approximate  formula 

a  [ - - - - - -i - - - - - £-]16  (b) 

c  +  a.  x  +  a_  x  +a„  x  +a.  x  +a_  x  +a,  x 
1  2  3  4  5  6 

where  x  =  - - 

7.  0710678119 

1 

c  =  1.0442737826  =  ifc, _ 

\TT 

ax  =  .03682270091 

a  =  .  11038499229 

a3  =  .  12101231517 

a  .  =  .  00992153425 
4 

a5  =  .  09025371127 
a6  =  .  07026624580 

for  (x)  in  the  range:  -  7.  0710678119  £  x  <  7.  0710678119 
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The  formula  (b)  will  yield  6-place  accuracy.  Differentiating  (a) 
we  obtain: 


P'(x) 


Now, - 1=  =  .  398994  22811 

*s/I7 


(c) 


Let  x2/2  =  x 
then 

— x  -  -2  -3  -4  -5  -6  -10 

e  =  a  -a  x  +  ax -ax  +a.x  -  a,  x  +  a,x  ....+alftx 
0  1  2  3  4  5  6  10 


(d) 


where 


a0=1*° 

a1=-1.0 

a2  =  .5 

ag  = .  166666 

66667 

a  = .  416666 

4 

66667 

a5  =  .  008333 

33333 

a6  =  .  001388 

888889 

a?  =.000198 

41270 

ag  =.000024 

80159 

a9  =.0000027 

5573 

a1Q=. 0000002 

7557 

(d)  will  yield  10-place  accuracy  if  the  following  identity  is  used: 


e  X  =  (e  n  )n  ,  where  ^  <  .  5  for  any  x  >0 
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(For  example,  if  it  is  desired  to  compute  e  ,  first  compute  e  ^  , 
then  raise  to  the  fourth  power). 

Now  we  can  compute  P(x)  and  P'(x) 

Let  f  (x)  =  P(x)  -  a 
f'(x)  =  P'(x) 

Iteration  formula  for  x: 

x  =*  .  .  x  .  (el 

n+1  n  f'(x)  n  P'(x) 

Take  an  arbitrary  value  for  x.  We  calculate  P(x)and  P'(x)  using 
(b),  (c).  Substituting  these  into  (e)  give  us  a  second  value  for  x. 

Repeat  iteration  by  repeated  application  of  (e).  The  process  converges 
to  the  desired  value  of  x.  If  this  value  is  y,  we  set: 

y  =  (  r-1)  sfiTj’ 

or,  r  =  (  y/"/n1)+  1 

Then,  for  a  given  n^,  r  is  computed 

II.  Computational  procedure  for  finding  r1,  given  a  ,  an  assigned 
value  of  P  where  P  is  given  by  the  formula  on  page  44.  As  before, 

S  8 

set  P  =  a.  The  incomplete  beta-function  I  (p,q)  can  be  represented 

8  X 

by  the  following  series: 

Ix(5.  5)  =  x5[a4x4-a3x3  +  a2x2  -a^  +  a^  (f) 


78 


where  a  =  70 
4 

a  =  315 
a2  =  540 
a^  =  420 

ao  *  126 

This  function  also  represents  P  from  formula  given  on  Page  45. 

s 

Now  f  (x)  =  a  -  1^(5,  5) 

f'(x)  =  -630  x4(l-x)4 

To  find  x,  we  use  the  following  iteration  formula: 


x  +  1  =  x 
n  n 


f(x) 

f'(x) 


=  X  + 

n 


a  -  Ix(5.5) 

4  4~ 

630  x  (1-x  ) 


Using  the  solution  of  this  iteration,call  it  y,  we  set : 

2 

u 

y  =  ^*°e  ’  2  du 

where  t  =  (r'-l)'VN^ 

and  the  problem  now  reduces  to  solving  for  t.  But  this  is  exactly 
the  same  type  of  problem  as  in  (I  ),  Using  the  method  given  in  (I) 
we  compute  t.  Then  for  a  given  we  compute  r'  from  the 
relationship  u  =  (r1- 1)  . 

III.  Computational  procedure  for  finding  Q  given  and  r. 

r  was  obtained  using  Computational  Procedure  I  from  page  46.  Q  is 
given  by: 
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n  ~  1 
Q  ~ 


du 


where 


To  compute  Q,  assign  a  value  to  Yj  (saY>  Yj  =  1*  1)  and  to  Nj  . 
Using  the  value  of  r  obtained,  we  calculate  u^  .  With  this  value  of 
Uj  we  use  formula  (b)  and  calculate  P(x) ,  where 


7.0710678119 

This  value  of  P(x)  is  the  desired  Q. 

IV.  Computational  Procedure  for  finding  Qgl  given  Y^»N^,r'.  The 

value  of  r*  was  obtained  using  Computation  Procedure  II.  Now  from  page  47 

2 

i  u2  .  JL_ 

9  /  e  2  du 

S  vZtt  -00 


where 


(r1  -  yi)^*1 


Assign  a  value  to  Yj  (say  Yj  =  1*  1)  and.  to  N^.  Using  the  value  of 
r’  obtained,  we  calculate  u^.  With  this  value  of  u^  we  use  formula 
(b)  and  calculate  P(x)  where  x=  7.  o7To^78119 

This  value  of  P(x)  is  q  .  In  formula  (f)  set  x  equal  to  the  value  of 

s 

q  which  we  have  obtained,  and  calculate  I  (5,5).  This  value  is  the 
s  x 

desired  value  of  Q  . 

s 
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V.  Given  the  function: 


${x)  = 


k/Ttt  -X 


/X  e  2  du 


4  (x)  can  be  approximated  to  six  decimal  places  by  the  following 
series  in  the  range  -  7.  0710678119  <  x  <  7.0710678119 


*<x)  =  1  -  [ 


1 


.16 


-2  -3  .4  .5  .6 

a„  +  a,  x  +  a„  x  +  a_  x  +a.  x  +a_  x  +  a,x 
0  1  2  3  4  5  6 


where  IT  7.  0710678119 
a0  =  1 

at  =.  03526  15392 

a  =  .  105705  03075 

a3  =  .  115881  7900 

a.  =  . 009500  89375 
4 

a5  =.  086427  2500 
a6  =  .067287  1875 
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FUNCT 
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f 


GRAPH  OF  PROBABILITY  OF  ERROR  AS  A  FUNCTION  OF  N 
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GRAPH  OF  PROBABILITY  OF  ERROR  AS  A  FUNCTION  OF  N 
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Fig.  7  GRAPH  OF  PROBABILITY  OF  ERROR  AS  A  FUNCTION  OF  N 
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ig.  9  GRAPH  OF  PROBABILITY  OF  ERROR  AS  A  FUNCTION  OF  N 
FOR  8  =  0.  9  ,  v,=  1.  5 
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GRAPH  OF  PROBABILITY  OF  ERROR  AS  A  FUNCTION  N, 


MINI 


LOWER  95%  CONFIDENCE  LIMIT  FOR 


n  (OBSERVED  NUMBER  OF  GRAINS  IN  TEST  AREA) 


Fig.  14  LOWER  95%  CONFIDENCE  LIMIT  FOR  N}  (  Vj  >  1 ) 
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MINIMUM  EXPECTEO  NUMBER  OF  GRAINS  IN  TEST  AREA  TC  SATISFY 


.2  .4  .6  .8  1.0  1.2  1.4  16  1.8  2.0  2.2  2.4  2.6  2.8  3.0 


Fig.  15  MINIMUM  EXPECTED  NUMBER  OF  GRAINS  REQUIRED  IN 
TEST  AREA  TO  SATISFY  P=.05,  Q  =  .05  FOR  /  1 
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n 


Fig.  16  95%  CONFIDENCE  INTERVALS  FOR  N  (v^1) 
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X 
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I- 


Fig.  18  NORMAL  APPROXIMATION  TO  THE  POISSON  DISTRIBUTION 
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Fig.  20  NORMAL,  APPROXIMATION  TO  THE  POISSON  DISTRIBUTION 


> 


Fig.  21  NORMAL  APPROXIMATION  TO  THE  POISSON  DISTRIBUTION 


Fig.  22  NORMAL  APPROXIMATION  TO  THE  POISSON  DISTRIBUTION 
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Fig.  2  3  NORMAL  APPROXIMATION  TO  THE  POISSON  DISTRIBUTION 


i 


i 

fcO 

u 


B 

d 

p. 

cd 

b 

o  b 

P  P 
O 


P.  B 
co  © 

O  P 
P  to 

o  a 


*4 


>» 

© 


o  a  o  g  d  a 

(gf  s-S  3  3 


3vT> 

ON  * 

i  Shh 

p  4  p*C"- 

^  6  3 

5  a  g  £  $ 

\°x£* 

jj  o3  J 

p  p  :3  gw 
b  d  r*  d 
fc.  O  J  •' 
o  J  b  i 

•  2  •  o  l 


a 

iT\  G 

-  ~  •  oND 

J©  3  a  p  oi 

B  •  d  I 

d  >«  a  a  o 
a  J  po 
O  d 

o  c  •  « 

O  O  CS  B  F-« 

.85  ^ 

M  u  U  w 

H  -H  Q  y  • 

C  -4  TJ 

o  a  3  o  r- 

•H  BOIS  S  O 

+>  o  m  cr 

3  a  h  •  j 

Sn  J  3  •  >• 
«d  •  a 

O  P  J  ^  d 

p  p  -4  dw 

©  JB  >  >.0  0 

•4  ft»  ^  JO  P  fP 


i  d 
H  P 
O  O 

•  a a 


44  “ 

BOB 

—  P  O 

O  d  o  >.«*' 


d  d  o  d  o  >»•*■  o) 

ar 44 


_  as 

P  •  ©  _ 

c  «  £  d 


, .  j  p  p  6 

d  p  p  P 

•a§  §  3  ^ 

O  Oi  O  B 

O  d  d 


•  P  P  ■ 

_  o  £  ■§)'§)  © 

g  HBSS  8 


d  d 
B 


H  d  S  p  3 

GPP  © 

o  p d 
o  <p  G  G 

O  O  d)ri 

b  d  b 

©  B  B  G  O 

d  ©  G  ©  o 
a  k  r<  o.h 
p  3  ci  o 

3*S^|5 

BO  -3  P 
P  O  © 

«  fl  S 

©  Otp  rl 


P 

•  G  •  < 

_  w  a  d  (Oi 

_  _  opqpud  •  t 

5  H  p  odu^© 

OB  d  ti  <0  I 

B  O  B  P  0*1 
©P©P.©©PG>»i 
dBddh^w-^p 
O  B  O  B  ©  3  P  < 

p  Cl  tf)  ©  •>  g  to  1 

©  -»  d  p  .  ©dr 

P©pp«p©^© 
35  ClTj  fl  d  ©  d 

U  P  fl  P  <H  P  u 
©  ©  ©  o  Pi 

Od  OQ.TiGGd' 
opo©c:(3paf3 
(4  <d  p  3  op 

a  m  -d  B  -d  M  A 

G  5  p  G  q  © 

CPdOP©dG© 

o  c  p  p  d  d  o  od- 

P  P  Q  O  P  P  P  p  ' 

p  o  3  «  p  © 

B  V.  ?  B  B  P  O 

3  P  ©  ©  G  3  p 


P  *  B 

o  d  P  G  S  p  *5 
P  d  O  P  B  ©  »d 
P  „  *0  p  a  p  o  a 
©  ra  p  p  B  *.a  o  p 
~  G  ©  P  «  B 


d  3  d  o  ©  d  © 
p  d  P  B  Di  d  P  B  (4 


#  ^  2 

p  3  o  o 
o  d  p 
b  p 

i  B  B  B  © 

»  a  p  u  © 

•  B  d  o  o  o 
i  ©  ?2>  u  o  p 

I  P  P  fi  M 

©  ©  Q.  © 
i  d  s  2 
-  p  v*  ©  «d 
d  o  d  • 
.  «  ©  p  ©  d 

l  O  B  B  O  © 

>  c  q  ts  d  d  n 
i  b  d  d  p  B  p 

i  P  O  p  p  H 

■  a  s « 

|  S  >J1 

: » 
i  p  d  ©  n  h  d 
I  O  H  P 
3  -d  O  p 

>  d  >.  B  •  B 

)  B  P  ©  B  d  3 

«  p  k  ©  o  a 
>«  ©  b  ©  d  p 
3  g  d  a  p  b 

ip  ©  H  o  B  73 

3  B  d  -d  ti  fl 

; >  *  .  s  s  a 

j#So  M  S)  d 

4  p  p  P  ©  P 

:  o  p  p  d  G  o 

>  ©  q  p  p  d 
0*  a  G 

s  x  d  3  o 

>  ©  p  <h  P 


B 


,  '*8 

.  >.  >>  3 

P  P  O 

P  P  f-t 

P  ©  B  B  U 

k  ft  ti  ti  m  3 

3  o  ©  ©  b 

o  o  xi  'd 
op  ^ 

B  d  •  p  P  • 

^  ©  d  d  B 

©  B  (4  iD  '*0  © 

a  ti d " ,H 


d  B 
G  o 
^  v  O  o 

ra  k-H  ftH 

p  3  en  o 

a'g 

BO  P  P 
P  0  ©  _ 

«  kd  fl  ^ 
P\  Oi  p  p  d 


&B  P  ©  * 

©  ©  p  d  >.  1 
‘  *  p  p 
3  p 

M  (0 

©  G  « 
P  ©  p  © 
p  ©  d  ©  d 
w  p  ft  p  U 
©  o  p 
Gi  ©  g  G  d 
©  d  9  h  d  w> 
d  P  3  op 

B  TJ  P  - 

U  G  G  © 

◦  P  ©  3  G  © 

p  d  d  o  o  d  • 

O  P  P  tfti  p  f 

BP  © 

ft  ?  a  to  k  o 

SG  3  P 
p  d 

^5-d  . 

H  P  H  ’-J  1 

©  P  ©  k  ® 

5  W  P  O  ftP 

o  ©  d  © 
d  Pkd  p  «  P 


•  p  . . 

>»  .  d  © 

p  >»  >  d 

!  P  P  P 

IP© 

I  p  d  d 

i  ©  p  p  a  © 

d  d  b  K 

i  b  d  © 

■  p  d  G  >» 

O  BP  © 

)  U  HP 

P  Pi  ©  W  © 

IP  3  G  d 

d  p  ©  fH 

i  ©  G  B  d 

•  d  p  ► 

p  p  • 

I  «h  B  ©  O 

i  0  P  P  P  P  © 

o  p  p  d 


©  q  p  p  d 

Pi  ©  q  o 


£r —  —  —  -  — - — - -iir  —  —  —  - 


I 

I 


1 


